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ABSTRACT 


This paper presents the results of optical studies of the crystals of long chain normal 
paraffins. The studies were restricted to crystals grown from solvent solution at ordinary 
temperatures; polymorphic crystal forms obtained under other conditions are not con- 
sidered. 

It is shown that all long chain -paraffin crystals, including very pure synthetic com- 
pounds and those occurring in nature, are optically active and invariably dextrorotatory. 
It is further shown that the crystals are normally orthorhombic, approaching hexagonal 
symmetry. No monoclinic crystals have been found. 

The cause of the optical activity is discussed. The writer holds that the zig-zag type of 
molecular chain, currently believed to exist in the long chain n-parafhins, cannot cause the 
observed optical activity. A spiral molecular chain for these compounds is proposed as a 
working theory. Such a spiral chain is believed to account for the observed optical rotation. 


INTRODUCTORY STATEMENT 


The present paper deals with the crystal] structure of normal paraffins 
within the wax range: that is, the homologous series of ”-paraffins above 
CisHge, which are waxy solids at ordinary temperatures. The observa- 
tions here summarized, and which are the basis of the paper, were 
made on -paraffin crystals ranging from C2; to C32, inclusive. These 
were crystallized on glass microslides from a solvent solution at ordinary 
temperatures, between 20°C. and 25°C. The crystals studied ranged 
from very pure individual w-paraffins to mixtures of two or three m- 
paraffins differing by one or a few carbon atoms. All materials studied 
were “‘pure”’ to the extent that non-paraffin impurities were absent, or 
present in such small quantities as not to affect materially the crystalli- 
zation. 

All optical observations were made with a petrographic microscope, 
and except where stated otherwise, were made at ordinary air tempera- 
tures. The light source was a 10 watt lamp with Corning “Daylight” 
filter. Observations and measurements of crystals using controlled 
higher-than-ordinary temperatures were made with a Leitz-Weygand 
heating stage, the temperatures being estimated to the nearest 0.1°C. 
with a probable error in most observations of +0.2°C. Crystal interfacial 
angles were determined by averaging at least five readings, the probable 
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error of a single reading being about +0.12°, and of the mean, about 
+.0.05°%, 

The primary objective of this paper is to present the results of the 
writer’s microscopic observations of #-paraffin crystals and his inter- 
pretation of the molecular and crystalline structures of the long chain 
normal paraffin series, based on these observations. These studies have 
been carried out intermittently over a period of four years, and have 
involved the preparation and examination of more than two thousand 
microslides of paraffin waxes from ninety five different sources. It is the 
writer’s belief that the results of these studies are worth recording, if 
only to attempt to clear up some of the apparently contradictory obser- 
vations to be found in an extensive literature on the subject. Further- 
more, certain optical properties of -paraffins described in the present 
paper are difficult if not impossible to explain if the molecule of the long 
chain -paraffin has the generally accepted zig-zag form, as interpreted 
from x-ray studies. 

To prevent possible misunderstanding, it is emphasized that the n- 
paraffin crystals described in this paper are exclusively those grown ona 
glass microslide from a solvent solution (carbon tetrachloride, benzene, 
toluene, or xylene) at ordinary temperatures well below the melting 
points. The same compounds, crystallized under different conditions, 
show quite different crystalline forms. This fact is largely responsible for 
the confusion one finds in reviewing the many published descriptions of 
paraffin crystals. The present paper does not attempt to cover the sub- 
ject of polymorphism in v-paraffin crystals, or the causes thereof. 


EXTERNAL CRYSTAL FORM OF NORMAL PARAFFIN 


The form and optical orientation of the typical x-paraffin crystal are 
shown in Fig. 1. The crystallization is generally orthorhombic, but ranges 
from orthorhombic to hexagonal. In the common form the pinacoid 
{001} and the rhombic prism {110} predominate. The pinacoid {010} 
is occasionally developed. The pinacoid {100} occurs rarely, but has not 
been observed where the pinacoid {010} is developed. 

The crystals are small. The long, or b-axis direction of the rhombic 
plate rarely is found to exceed 0.7 mm., and the usual length for the 
larger crystals on the average slide is from 0.1 mm. to 0.5 mm. The av- 
erage thickness, in the c-axis direction, is of the order of 0.005 to 0.01 mm. 

Euhedral crystals of a very pure! -paraffin are nearly perfect in form, 


* Normal paraffin wax compounds of 100% purity probably do not exist. Paraffins 
graded as “‘c.p.”” may contain up to several per cent of other members of the same homolo- 


gous series. The term “‘very pure’’ as used throughout this paper indicates a compound 
probably at least 98% pure. 
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the corresponding interfacial angles in all crystals being nearly the same. 
Mix crystals of two or more w-paraffins, or crystals of a single -paraffin 
containing impurities, are generally more or less distorted and show con- 


(110) 


Fic. 1. Form and Optical Orientation of the Orthorhombic n-Paraftin Crystal. 

In its simplest form, as illustrated, the crystal is a thin rhombic plate. In the figure, the 
thickness of the plate is exaggerated. It should be noted that since the angular measure- 
ments are made with the microscope and depend on the interior angles of the rhombic 
plate, all descriptive data are given in terms of these angles. The obtuse interior angle is 
designated as 6, and the acute interior angle as ¢. Thus @ is equivalent to the conventional 
interfacial angle (110) \(110), while ¢ is equivalent to (110) (110). 


siderable variation in interfacial angles, both within the same crystal and 


among different crystals in a single preparation. 
The acute interior angles (¢ in Fig. 1) are more variable within a single 
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crystal than are the obtuse interior angles, @. In crystals of very pure 
n-paraffins, ¢1=¢2 (approximately) and 0,=6 (approximately). In mix 
crystals and in crystals of impure n-paraffin, as a general observation 
o1% do, whereas 6,=6» (approximately). In other words, the b-axis of the 
crystal tends to remain an axis of symmetry; the a-axis does not. 

The range in variation of ¢ angles and the relative frequencies of oc- 
currence of different values of ¢ within the observed range are shown on 
the chart (Fig. 2). This chart shows the variations found from about 800 
individual measurements of ¢ angles made on approximately 400 differ- 
ent crystals. To interpret this chart correctly, it should be kept in mind 


5S 60° 65° 70° Go 


Fic. 2. Frequencies of Occurrence of Various Values of m-Paraffin Crystal Angles (¢). 
Data: 800 measured ¢ angles on 400 crystals from 10 different sources. The angles were 
measured to the nearest 0.1°+0.05°, and were plotted on this chart to the nearest 0.5°. 


that the relative frequencies of ¢ angles shown pertain only to the 400 
crystals from 10 different sources which were measured to obtain the data 
for this particular frequency curve. Similar measurements on crystals 
from other sources, or on a much larger selection of m-paraffins from all 
natural and synthetic sources would probably show frequency curves dif- 
fering from that in Fig. 2. 

It should also be kept in mind that the frequency curve in Fig. 2 shows 
the frequencies of individual ¢ angles; not the means of ¢; and ¢» for the 
different crystals. In a single crystal the difference between ¢; and ¢» 
has been found to be as great as 20°, though the average difference is 
much less. 

Several interesting features are evident from an inspection of Fig. 2. 
The frequency curve declines rapidly to 60° (the angle for a hexagonal 
crystal form), and remains low for all values less than 60°. Frequency 
peaks, in order of importance, are seen at 68°, 65°, 66.5°, 64°, 74°, and 
71°. The peak at 74° represents very pure Octacosane (CosHss). The re- 
mainder of the data incorporated in the chart are from mix crystals of 
refined n-paraffins from petroleum sources. 
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In general, the entire range of ¢ values will not be found in mix crystals 
from any single source. Each sample of paraffin from a different source 
shows a characteristic frequency curve of ¢ angles, usually distinguish- 
able in its form from $-frequency curves of samples from other sources. 

Numerous observations of both pure and mix crystals show that the 
range of ¢ values in a pure individual m-paraffin is quite limited, and 
much greater in the mix crystals. For example, a very pure 2-Octacosane 
showed a sharp ¢-frequency peak at about 74.25°, and a range of ¢ values 
from 73.5° to 75.5°. A sample of pure u-Dotriacontane (C3:H¢s) showed 
a peak at 68° and a range from 65.5° to 69.5°. An equi-molar mixture of 
pure Octacosane and fairly pure Nonacosane showed a broad peak at 
68° and a range from 63.5° to 75.5°. Typical samples of mix crystals, re- 
fined from petroleum sources, showed one or more frequency peaks, usu- 
ally not sharp, and ranges in ¢ values in some cases as great as 15° or 20°. 

The evidence indicates clearly that the spread or range in ¢ values for 
any given sample is related to the purity of the material, or to the num- 
ber of individual v-paraffins present in the mixture. Although no exact 
relationships have been determined, it may be stated that the spread of 
¢ values varies inversely as the degree of purity of the material crystal- 
lized. With carefully controlled molar mixtures of very pure ”-paraffins it 
may prove possible to establish a mathematical relationship between the 
interfacial angles of the crystal and the molar percentages of -paraffins 
comprising the mixture. The writer has been unable to make such a 
study because he has been unable to obtain samples of more than one 
n-paraffin of sufficient purity. 


TWINS 


The long chain -paraffins crystallized from solution at ordinary tem- 
peratures, exhibit a variety of twins. The most commonly observed 
twinning planes are {100}, {110}, {310}, and {010}, listed in order of 
frequency of occurrence. Twinning apparently occurs in other directions, 
but much less frequently. 

Micro-twinning is apparently responsible in many instances for curved 
{110} faces observed in -paraffin crystals. In the majority of these in- 
stances, the curvature of the {110} faces is convex outward, but the re- 
verse curvature is common in crystals of m-paraffin from certain sources. 

Curved faces may also be produced where impurities are present in 
the solution, or in mix crystals containing several m-paraffins. Curvature 
is also common in crystals of the lower melting point members of the n- 
paraffin series when the temperature during crystallization is within 8° 
or 10°C. below the melting point. 

It is sometimes difficult to determine whether the curved faces are the 
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result of twinning, or of warping of the crystal structure produced by 
impurities, or of other factors. It has been observed that freshly formed 
perfect crystals may in the course of time develop curved or irregular 
faces. This is apparently the result of resorption produced by relatively 
low melting point amorphous wax or oil, which makes up the residual 
matrix surrounding the paraffin crystals after the solvent has evaporated. 
This feature is particularly noticeable where relatively low melting point 
crystals have not been stored in a cool place. 

In addition to twinning, a common feature of -paraffin crystals is the 
formation, by parallel growth on the {001} face, of one or more rhombic 
plates of successively smaller dimensions. These are usually more or less 
concentrically arranged, so that the entire crystal has the appearance of 
a very flat step pyramid. 

Occasionally, and particularly in pure m-Octacosane, the successive 
steps or tiers of rhombic plates, instead of having a concentric arrange- 
ment, are superimposed in such a manner as to present the appearance 
of a spiral arrangement,’ and in such cases, the pseudo-spirals may ap- 
pear right-handed of left-handed, with no statistical evidence to indicate 
a tendency toward either direction. 


BEHAVIOR IN PLANE POLARIZED LIGHT 


In plane polarized light the pure crystals are colorless, translucent to 
transparent. They show no sign of pleochroism. When crystallized from 
solution on the microslide beneath a cover glass, the rhombic plates are 
so thin that a reduced or an oblique light is usually necessary to make 
them visible. The refractive indices of the crystals are commonly only 
slightly different from the index of the surrounding matrix. This matrix 
is composed of amorphous waxy materials, and shows few or no indica- 
tions of crystalline structure. 

In the case of extremely pure single members of the -paraffin series, 
however, practically the entire area of the microslide preparation is crys- 
talline. In such cases, it is unusual to observe euhedral crystals, formed 
free of interference from neighbors, and the true boundaries of crystals 
are difficult to observe without the use of crossed nicols. 


BEHAVIOR BETWEEN CROSSED NICOLS 


Between crossed nicols in parallel rays (without use of condenser lens) 
the crystals behave differently, depending upon whether the preparation 
is of a single very pure -paraffin on the one hand, or on the other a mix- 


* This form of parallel growth with a spiral pattern differs from the spiral structure in 
certain paraffin crystals described by Prof. Chas. M. Heck, North Carolina State College, 
Raleigh (Heck,1936). 
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ture of two or more pure m-paraffins, or from one to several n-paraffins 
containing impurities. In general, pure preparations of a single m-paraffin 
compound show fairly marked double refraction in the c-axis direction, 
whereas the mix crystals and the impure preparations are usually feebly 
anisotropic, or practically isotropic (uniaxial). 

Interference colors are invariably first order gray, even in the pure 
unmixed compounds. This is true for the thickest crystals that have been 
observed. The insertion of a gypsum plate effectively enhances the colors. 

The following characteristics in parallel light between crossed nicols 
have been observed without exception in all ”-paraffin crystals which are 
appreciably biaxial: 


1. The faster ray vibrates approximately parallel to db. 


2. Extinction is approximately symmetrical, but as will be described later, it is in some 
instances noticeably asymmetrical, and the extinction angle varies slightly with dif- 
ferent crystals. 


3. Double refraction characteristically varies within a single crystal where the thickness 
in the c-axis direction is variable. Many crystals are thicker in the center and thin out 
toward the margins, being roughly of the form of a flattened pyramid. In such crys- 
tals, the central portion may be rather strongly anisotropic, and the marginal zone 
apparently isotropic. This observation suggests that many crystals appear isotropic 
because of their very slight thickness. 

4. When the b-axis of the rhombic plate is rotated into the “northwest” quadrant (with 
reference to the vibration planes of the nicols) the crystals show stronger double 
refraction than when the b-axis is rotated to the “northeast” quadrant. In thick, 
strongly anisotropic crystals, or in crystals of a single very pure n-paraffin, this 
phenomenon is not conspicuous, and might be missed if one were not looking for it. 
In the more feebly anisotropic crystals, and particularly in very thin crystals, this 
feature is striking. Such crystals, during a complete rotation of the microscope 
stage, may show double refraction only twice; both when the b-axis of the crystal 
is in the “northwest”? quadrant. When the b-axis is in the “northeast”? quadrant, 
there is either a very faint appearance of double refraction, or the crystal remains 
dark. There is nearly always complete extinction when the b-axis is parallel to the 
plane of the lower nicol. When the b-axis is parallel to the plane of the upper nicol, 
extinction is often incomplete. When the b-axis is oriented with the plane of the lower 
nicol, extinction is almost always symmetrical. When the a-axis is oriented with the 
plane of the lower nicol, however, extinction is often observed to be slightly asym- 
metrical. Wherever the angular positions of extinction and of maximum double 
refraction have been measured carefully, they have been found to be shifted or ro- 
tated slightly in a counter-clockwise direction from their normal positions (with 
the possible exception, previously noted, when the b-axis is parallel to plane of lower 
nicol). 

In order to investigate more thoroughly the nature of this anomalous behavior, 
the microslides were in many cases placed up-side-down on the microscope stage, 
and the identical crystals observed again in the inverted position. The phenomenon 
remained unchanged, with maximum double refraction still in the “northwest” posi- 
tion, and the positions of extinction and of maximum double refraction shifted 
slightly counterclockwise. 
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Breas Fic. 4 


Fic. 3. Photomicrograph in plane polarized light of a typical n-paraffin mix crystal 
of about 56°C. melting point, occurring in petroleum (Panhandle, Texas). Magnification 
<415. 

Fic. 4. Photomicrograph in polarized light of crystals of a pure synthetic n-paraffin 
(n-Dotriacontane, C32 H¢s). Crossed nicols, with parallel light. Magnification 48x. Vi- 
bration planes of the nicols are parallel to edges of the photograph; P-P’ being ‘‘north- 
south” and A-A’ being ‘‘east-west.”’ . 

The crystals seen in the photograph are all oriented with their c-axes “‘vertical”’ (nor- 
mal to plane of the paper) and with their b-axes (long diagonal of the rhombic plate) ori- 
ented at random in the “‘horizontal”’ plane (plane of the paper). It will be observed that 
those crystals with b-axis aligned in the ‘‘northwest’’ quadrant show relatively strong 
double refraction, whereas those with b-axis parallel to vibration plane of either nicol or 
aligned in the “‘northeast’’ quadrant are dark, or show relatively weak double refraction. 
This anomalous behavior is explained as an effect of the optical rotatory power of the n- 
paraffins. 


BEHAVIOR IN CONVERGENT POLARIZED LIGHT 


The interference figure is biaxial positive. Because of the invariable 
orientation of the crystals on the microslide with the c-axis in the vertical 
position, the interference figure always seen is the acute bisectrix figure. 
The isogyres show no accompanying isochromatic curves, due to the rela- 
tively small double refraction (8-a), and to the thinness of the crystals. 

Mix or impure crystals show small optic angles, and in many cases 
the interference figure appears to be practically uniaxial. In pure un- 
mixed n-paraffins, values of 2E have been observed as large as 41°. The 
usual range of 2E for distinctly biaxial crystals is between 20° and 40° 
Among the mix and impure crystals, the 2E angles range generally from 
0° to 10° or 15°. Both in the strongly biaxial and in the approximately 
or apparently uniaxial crystals, the sign of the interference figure is in- 
variably positive. No n-paraffins of opposite sign have been observed. 

The writer has been unable to determine accurately the refractive in- 
dices of -paraffin crystals in the three crystal directions. In the in- 
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variable orientation of the crystals on the microslide, only the 8 and a 
indices may be measured conveniently, and these two indices are usually 
nearly equal. Measurements of these have shown a range of values from 
1.50 to 1.52. If the average of @ for all the observed paraffins is assumed to 
be 1.51, then the observed range for distinctly biaxial crystals of 2E = 20° 
to 40° would correspond to the values 2V = 13° to 26°. 

The interference figure is always central in the field of view, and shows 
no displacement when the microscope stage is rotated. Thus the acute 
bisectrix is parallel to the optical axis of the microscope, and normal to 
the (001) face of the crystal. No evidence of the slightest inclination of 
the acute bisectrix has been observed in any of the many n-paraffin crys- 
tals from all sources, both natural and synthetic. Because of the impossi- 
bility of thus detecting a slight inclination of the acute bisectrix, these 
observations do not necessarily prove that every crystal studied was 
either orthorhombic or hexagonal. The observations do, however, show 
that none of the strongly monoclinic crystal forms reported from x-ray 
studies, are to be found among the long chain n-paraffins crystallized 
from solution at temperatures well below transition and melting points. 

As would be expected, in view of the previously described asymmetri- 
cal behavior of -paraffin crystals between crossed nicols with parallel 
light, the behavior in convergent polarized light is also anomalous. 
Without any observed exception, even in the purest -paraffins, the 
axial angle (2E) is definitely larger when the b-axis of the crystal is 
oriented in the ‘‘northwest”’ quadrant. This asymmetry is usually but 
not invariably more pronounced in those crystals showing 2E values 
somewhat below the average for all crystals studied. The degree of the 
asymmetry is fairly constant for all crystals from any one sample or 
source. This statement is particularly true of pure, unmixed m-paraffins. 

Careful measurements of the axial angles of 24 crystals from 15 differ- 
ent sources, ranging from very pure -paraffins to mixtures and impure 
preparations, and including both natural and synthetic n-paraffins, 
showed the following results: 

2E(N.W.)—Range 17.8° to 41.6°; arithmetic mean =30.2°. 


2E(N.E.) —Range 15.3° to 38.9°; arithmetic mean=25.0°. 
DEN Wear NE) oct tsteteraroere arithmetic mean =27.6°. 


Maximum asymmetry measured in a mix crystal; 2E(N.W.) =30.9°; 2E(N.E.) =20.4°. 
Minimum asymmetry measured in a pure crystal; 2E(N.W.) =40.7°; 2E(N.E.) =38.9°. 
Minimum asymmetry measured in a mix crystal; 2E(N.W.) =22.1°; 2E(N.E.) =20.4°. 


All determinations of 2E on individual crystals represent means of 
repeated measurements. The probable errors of the means are in most 
cases less than +0.5°. 
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INTERPRETATION OF ANOMALOUS BEHAVIOR IN POLARIZED LIGHT 


The previously described anomalous behavior of -paraffin crystals 
between crossed nicols shows that these crystals are optically active. The 
observations furthermore show that all the crystals thus far examined are 
dextrorotatory. 

The dextrorotatory action of the #-paraffin crystal is illustrated graph- 
ically in Fig. 5. The left diagram (Fig. 5a) represents a crystal with the 
b-axis oriented at 45° in the ‘‘northwest”’ quadrant; the right diagram 
(Fig. 5b) another crystal with the b-axis oriented at 45° in the “north- 
east” quadrant. In both diagrams lines PP’ and AA’ represent the vibra- 
tion directions of the lower and upper nicols respectively. In each dia- 


P 


Fic. 5. Optical Rotation of n-Paraffin Crystal Between Crossed Nicols. 


gram the b-axis of the crystal is shown by a solid line; the a-axis by a 
dashed line (cf. Fig. 1). The projection of the vertical c-axis on each 
diagram at the intersection of the a and 0 axes, and is designated by the 
letter C. The c-axis of the crystal on each diagram coincides with the 
optical axis of the microscope. The distances CB =CB’ represent the 
amplitude of the plane polarized ray emerging from the lower nicol. 

When this polarized ray enters the crystal, it is split into two polar- 
ized rays, the ordinary, or (0) ray and the extraordinary, or (e) ray. If the 
n-paraffin crystal were optically inactive, the (0) ray would vibrate 
through the crystal in a direction parallel to the b-axis, and the (e) ray 
would vibrate parallel to the a-axis. 

An optically active dextrorotatory crystal would produce a clockwise 
rotation of the vibration direction of one or of both rays. Observations of 
n-paraffin crystals indicate that the amount of rotation suffered by the 
(0) ray is usually small or negligible. For greater simplicity, therefore, 
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the (0) ray in Figs. 5a and 5d is assumed to have undergone no rotation. 
The clockwise rotation of the (e) ray is shown by a solid line in both 
diagrams. The amount of rotation is about 25° in Fig. 5a and about 20° 
in Fig. 5d. 

At the 45° position the component amplitude of the (0) ray which 
emerges from the crystal is represented by the distance C-o, found by 
dropping a perpendicular from B to the b-axis. A perpendicular to AA’ 
is then dropped from the point (0) to point (oe’) on AA’. Then the 
distance C-oe’ represents the component amplitude of the (0) ray which 
emerges from the upper nicol. 

The (e) ray, which like the (0) ray, originates the instant the plane 
polarized light enters the crystal from below, vibrates initially parallel 
to the a-axis (dashed lines in Figs. 5a and 50). The component amplitude 
of the (e) ray while passing through the crystal is represented by the 
distance C-e, found by dropping a perpendicular to the a-axis from the 
point B’. In passing through the crystal the plane of vibration of the (e) 
ray is rotated in a clockwise direction, so that on emerging from the 
upper surface of the crystal the ray has changed its vibration direction 
from C-e to C-e’. A perpendicular to AA’ dropped from (e’) to (ee’) gives 
the distance C-ee’, which is the component amplitude of the rotated 
(e) ray which emerges from the upper nicol. The sum of the components 
C-oe’ and C-ee’ gives the total amplitude of the combined rays emerging 
from the upper nicol. 

It may be visualized from Figs. 5a and 5d that as the dextrorotatory 
crystal is rotated on the microscope stage between crossed nicols, the 
sum of the components of the two emergent rays (oe’ and ee’) will vary 
greatly, and the pattern of the variation will differ from that of an op- 
tically inactive crystal. The pattern will also depend on the amount and 
direction of rotation of either or both rays emerging from the crystal. 

To illustrate the effect of optical activity, the two curves in Fig. 6 
show the changes in relative intensity of transmitted light as the crystal 
is rotated (between crossed nicols) through 180° in a clockwise direction, 
starting (on the left) with the b-axis parallel to the vibration plane of the 
upper nicol. The dashed line curve is that of an optically inactive crystal. 
The solid line curve is that of a -paraffin crystal in which the (0) ray 
has undergone no rotation, and the (e) ray has undergone a clockwise 
rotation of 20°. The curves need not be extended as they are repeated 
over each 180° rotation of the crystal. 

The solid line curve (for an optically active crystal) in Fig. 6 illustrates 
qualitatively several observations of the anomalous behavior of n- 
paraffins in parallel polarized light which have been described previously. 
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INTEHSITY OF LIGHT 


90° 45° 9° 45° 90° 
NORTH WEST QUADRANT NORTH EAST QUADRANT 
Fic. 6. Orientation of b-Axis of Crystal vs. Light Intensity. 
Dashed Line Curve; r=0° 
Solid Line Curve; r(o) =0° 
r(e)=+20° 
These are: 


1. Light intensity is invariably greater with the b-axis at 45° in the ‘‘northwest”’ quad- 
rant than at 45° in the ‘‘northeast’’ quadrant. 

2. Extinction is often incomplete and slightly asymmetrical with the b-axis in the 
“east-west’”’ position. 

3. Extinction is usually complete and symmetrical when the b-axis is in the “north- 
south”’ position. 

4. Positions of maximum light intensity and positions of extinction (where measurably 
asymmetrical) are always shifted counter-clockwise from their respective normal 
positions. 

5. In the case of a thin crystal, or one with weak double refraction, the light intensity 
maximum in the ‘“‘northeast’’ quadrant would not be observed, and the crystal would 
appear to remain dark through 180° of rotation on the microscope stage. 
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The solid line curve in Fig. 6 is intended to illustrate only qualitatively 
what actually happens. It is not intended to imply that only the ex- 
traordinary ray is rotated; and it may be assumed that both the (e) 
and (0) rays undergo amounts of rotation which vary with the structure 
and thickness of each individual crystal. However, observations suggest 
that the (e) ray is rotated more than is the (0) ray in the great majority 
of -paraffin crystals. 

In order to check further the theoretical asymmetry shown in Fig. 6 
against actual behavior of the -paraffin crystal in polarized light, a 
single crystal was selected, and a number of repeated measurements 
made and averaged. The crystal selected was a fairly large distinctly bi- 
axial mix crystal, the melting point of which was later found to be about 
61° C. The axial angles (2E) were 27.4° (N.W.) and 24.8 (N.E.). The 
results of the measurements are recorded in Table 1. 


TABLE 1 
; No. of | Prob- 
Obs. ; Approx. Revs nN oN repeated] able 
No. Light filter used length of incident (N.W.)] (N.E.) Sum eenturellerranor 
light-peaked at: eur Kaean 
1(a) | Red shade yellow | Above 610 mp | 48.0° 10 +0.6° 
Corning No. 3480 
1(b) | Red shade yellow | Above 610 mz 36.2° | 84.2° 10 +0.6° 
Corning No. 3480 
2(a) | Nene Approx. 589 my | 48.4° 10 020" 
2(b) | None Approx. 589 mz 43 .8° | 92.2° 10 a OFSe 
3(a) | Dark blue green | Below 520 mz 48 .3° 10 +0.7° 
Corning No. 4303 
3(b) | Dark blue green | Below 520 mp 45.0° | 93.3° 10 +0.4° 
Corning No. 4303 


Light source was 10 watt tungsten lamp with Corning “Daylight” glass. 


In this table, the angle \ is the angle between the plane of the lower 
nicol and the b-axis of the crystal when the crystal is rotated to the 
position of maximum light intensity in the quadrant indicated. 

The results shown in Table 1 agree qualitatively with the theoretical 
requirements illustrated in Fig. 6. Thus it is evident that in dextrorota- 
tory crystals the peaks of light intensity are displaced in a counter- 
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clockwise direction. Still another feature brought out in Table 1 is the 
indication that the longer the wave length of the transmitted light, the 
greater the counter-clockwise shift of the light intensity maximum in the 
“northeast” quadrant. 

It is obvious from the preceding observations that it would be possible 
to determine with fair accuracy the amount of optical rotation in a - 
paraffin crystal by measuring accurately, with monochromatic light, the 
positions of extinction and of light intensity maxima, and the relative 
peak intensities. To determine the specific rotation, the thickness of the 
crystal would have to be measured. Thus it would be possible, and of 
interest, to compare the specific optical rotation of a paraffinic petroleum 
or petroleum fraction with the calculated specific rotation of »-paraffin 
crystals extracted from the same petroleum. 


CAUSE OF OPTICAL ACTIVITY OF N-PARAFFIN 


It has been pointed out that whereas the observed optical asymmetry 
of the ”-paraffin crystal between crossed nicols is often quite pronounced 
in the impure and mix crystals, it is still clearly present in the purest 
synthetic m-paraffin crystals. It is therefore highly probable that this 
optical activity is caused by the internal structure of ”-paraffin, and not 
necessarily by the presence as impurities of some foreign organic com- 
pounds, as seems generally to have been assumed in the case of paraffinic 
petroleums. 

In the above connection, Gruse and Stevens (1942) state: ‘Optical 
rotatory power is exhibited by nearly all crude oils. The majority of 
them are dextro-rotatory; some are laevo-rotatory; and only a few appear 
to be inactive. The optical activity is not uniformly distributed through- 
out the distillation range of the crude but is usually a maximum for the 
intermediate fractions. Curiously, the maximum occurs at about the 
same boiling point for all crude types, 250° to 300° C. at 12 mm. Hg., but 
no extensive data have been collected on this point. No activity has been 
observed in fractions boiling below 200° C. at atmospheric pressure. The 
decrease in activity above the maximum may be partly due to racemiza- 
tion by heat. The existence of optical rotation of oils has given much 
weight to the organic theories of the origin of petroleum.” 

Of course a considerable number of organic compounds, known or be- 
lieved to be present in petroleum, could cause the observed optical ac- 
tivity. It is perhaps significant, however, that the range of crude oil 
fractions showing the maximum optical activity corresponds in a general 
way with the range of fractions which in paraffinic crudes contain most 
of the long chain m-paraffin waxes. The fact that the majority of crude 
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oils are dextrorotatory is likewise significant. Absence of optical activity 
in a crude oil might be due to racemization, or it might be due to absence 
of optically active constituents, such as long chain m-paraffins. The laevo- 
rotatory activity of some crude oils may be caused by the presence of 
laevorotatory hydrocarbons. At least, the evidence presented in the 
present paper would seem to show that the laevorotatory activity, ex- 
hibited by some crude oils, cannot be due to the presence of n-paraffins. 

It has been shown by Miiller (1928, 1929) that the molecules of the 
long chain v-paraffins are arranged with their long axes parallel to the 
c-axis of the crystal, and the ends of the molecules lie in planes normal 
to the c-axis. Miiller (¢bid.) and others (Piper et al., 1931; Kolvoort, 
1938; Gray, 1943) have also shown that under certain conditions at least 
some members of the -paraffin series may have the long axes of the 
molecules inclined to the c-axis and to the {001} planes. 

In the present paper it has been shown that when a-paraffin crystals 
form from solvent solutions at ordinary temperatures well below transi- 
tion and melting points, the molecular chains must be vertical, or parallel 
to the c-axis of the crystal, and that the crystals are primarily orthorhom- 
bic. Under the conditions specified, no definitely monoclinic crystals 
have been observed. 

It appears that nothing in the packing arrangement or attitude of the 
n-paraffin chains in forming the crystal unit cells and the lattice can ex- 
plain the observed optical activity. If the optical activity is not caused 
by the presence of foreign material in the inter-molecular spaces of the 
lattice, then it can only be produced by the internal structure of the 
n-paraffiin molecule. 

Miiller (1928) working with one or two perfect crystals of very pure 
n-Nonacosane, carried out the only complete x-ray study that has been 
made of any long chain m-paraffin. This work has become classical, and 
Miiller’s conclusions on the structure of the long chain -paraffin mole- 
cule and crystal lattice have become generally accepted. 

Briefly stated, Miiller found the m-paraffin molecular chain to be a zig- 
zag form, in which the carbon atom centers all lie in one plane, and in 
that plane, form two vertical rows on opposite sides of, and equidistant 
from the zig-zag chain axis. Instead of the carbon-carbon bonds forming, 
one with the next, the tetrahedral angle of 109° 28’, as might be expected 
from the known tetrahedral form of the carbon atomic domain, the zig- 
zag angle was found to be between 76° and 92°. The chain axis was found 
to be parallel to the c-axis of the crystal, and the zig-zag planes of the 
molecules were found to lie at an angle to the b-axis of the crystal of some- 
where between 23° and 30°. In any single {001} layer of molecules, half 
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the molecules have their zig-zag planes rotated 23°-30° clockwise; the 
other half counter-clockwise. In each {001} layer of molecules, the mole- 
cules are staggered along rows which are parallel to the a-axis. Each row 
alternately has molecules with zig-zag planes rotated clockwise or with 
planes all rotated counter-clockwise. 

The above described zig-zag molecule and crystal lattice (for ortho- 
rhombic crystal forms) have come to be generally accepted as represent- 
ing the true structure of the long chain -paraffin series. There is nothing 
in this molecular structure or in this lattice structure which can explain 
the observed dextrorotatory optical activity of m-paraffins, and it is 
therefore apparent that in some respects, Miiller’s reconstruction from 
x-ray data must be in error. In general, the results of the optical studies 
described in the present paper do not disagree with Miiller’s reconstruc- 
tion of the crystal lattice. Exception is taken only to Miiller’s reconstruc- 
tion of the molecule, on the grounds that the optical observations tend 
to confirm the essential accuracy of Miiller’s reconstruction of the crystal 
lattice, and that neither this lattice arrangement nor the zig-zag form 
molecule can account for the observed optical rotation. 


SPIRAL STRUCTURE OF THE LONG CHAIN N-PARAFFIN MOLECULE 


As a working theory, subject to confirmation by further investigations, 
it is suggested that the long chain v-paraffin molecule has a spiral struc- 
ture. 

Starting with this basic assumption, which is not new as applied to 
some other long chain compounds, certain other reasonable assumptions 
will be made. These are: 


1. The diameter of the carbon atom domain, or length of the C-C bond, is 1.54 A, as in 
diamond. 

2. Adjacent C-C bonds form the tetrahedral angle, 109° 28’. From this assumption, a 
spiral or solenoidal arrangement of carbon centers is a logical consequence. 

3. Of many different geometrically possible tetrahedral spirals, the two simplest ar- 
rangements are selected, and from these two, the left handed spiral is chosen. (See 
photo of model, Fig. 7.) 

4. The C-H bonds also form the tetrahedral angle with all other bonds from a common 
center. This assumption may be unnecessary to explain the optical rotation. It is 
made merely to complete the picture in an orderly manner. 


From the foregoing assumptions, the following relationships and di- 
mensions of the spiral chain (see Fig. 7) have been calculated: 


1. The carbon atom centers are spirally arranged around the chain axis, being equally 
spaced on a left-handed helix, the parametric equations of which are; x=R cos 6; 
y=R sin 6; z=K@, where R=0.534 A, K =0.5189, and 6=131° 43.94’ (=2.2991629 
radians). 
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2. Various dimensions and angles are: 
(a) Diameter of the helix=1.068 A. 
(b) Cross sectional area of helix=0.895 (A)2. 
(c) Inclination of C-C bonds to direction of chain axis=39° 14.035’. 
(d) Spacing (in direction of chain axis) of successive carbon centers=1.193 A. 


Fic. 7. Model of a spiral molecule of n-paraffin (-Octacosane). The centers of the 
tetrahedrons correspond with the carbon atom centers, which fall on a helical curve. For 
greater convenience in constructing the model, the carbon-carbon bonds are made to inter- 
sect the tetrahedron face centers rather than the corners. The helical arrangement of the 
carbon atom centers is the same in either arrangement. Positions of the hydrogen atoms are 


shown by the black dots. 


3. The molecule has no plane of symmetry. 

4. As is evident from the calculated value of 8, no two carbon centers in a chain of any 
reasonable finite length lie in vertical alignment. There is, however, an approximate 
repetition or cycle for every eleventh carbon atom. Thus, C-12 is situated almost 
vertically below C-1; C-13 below C-2; C-14 below C-3, etc. 
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A reconstruction of the assumed x-paraffin spiral molecule is shown in 
Fig. 7, which represents a chain of 28 carbon atoms (n-Octacosane). For 
convenience in construction, the C-C bonds are made to intersect the 
tetrahedron face centers, though they could as well be drawn through the 
tetrahedron corners. The result is the same as long as the bond length of 
1.54 A is retained. The hydrogen atoms are represented on the model 
by black dots at the centers of the ‘‘free’’ tetrahedron faces. It is seen that 


TABLE 2 


Measurement Miller Spiral Theory Remarks 


Length of C-C Bond | 1.8 to2.0 A | 1.54 A (exactly) 


Distance between 1.2 to1.6A | No rows exist Diam. of helix=1.068 A 
2 vertical rows : 


Dist. between con- | 2.537 A No rows exist Dist. parallel to chain axis be- 

secutive centers are tween alternate carbon centers 

in either row =2.386A 

Zig-zag angle 76° to 92° Tetrahedral angle | If C-C bondswere projected on 
= 109° 28’ plane parallel to chain axis, ap- 


parent angle would vary be- 
tween 109° 28’ and 180°. 


Angle between C-C | 44° to 52° 39° 14.035’ Projected on plane as per 
bonds and chain above, apparent angle would 
axis direction vary between 0° and 39° 14.035’ 
Angle between zig- | 23° to 30° non-existent 

zag plane and crys- 

tal b-axis 


in addition to the spiral arrangement of the carbon centers, there are six 
spirals made by the hydrogen atoms in the CH groups It is not neces- 
sary to assume that the hydrogen spirals exert any influence in producing 
the observed optical] activity. 

For comparison, Table 2 shows the spacings and angular relationships 
of the carbon centers as reconstructed by Miiller from «x-ray studies, 
and as result from the spiral arrangement described above. 

According to Miiller’s calculation of 2.537 A for the distance between 
carbon centers along either row of the zig-zag chain, the addition of one 
carbon atom would increase the chain length by 1.2685 A+0.5%. With 
a spiral chain of the dimensions described, the corresponding increase in 
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chain length would be 1.193 A. It should be noted, however, that Miil- 
ler’s value corresponds to a calculated C-C bond length of about 1.9 A 
whereas we have assumed for our spiral chain a bond length of 1.54 A. 
If we were to assume a bond length of 1.9 A for the spiral chain, the addi- 
tion of one carbon atom would increase the chain length by 1.47 A. 

Whether or not Miiller’s x-ray data could be re-interpreted to fit a 
spiral molecule is a question beyond the scope of the present paper. It is, 
of course, possible that the zig-zag type of chain exists in some m-paraf- 
fins, or in all m-paraffins under some conditions. On the other hand, the 
optical evidence presented appears definitely to favor the spiral form, as 
it is difficult to visualize how either the zig-zag chain or the lattice struc- 
ture could produce the observed optical rotation. Several thousand in- 
dividual »-paraffin crystals from many sources have been studied in 
polarized light. Among these not a single crystal was observed which 
failed to show the characteristic indications of optical activity. 


BEHAVIOR OF N-PARAFFIN CRYSTALS ON HEATING 
REVIEW OF PREVIOUS WORK AND THEORETICAL CONSIDERATIONS 


Members of the -paraffin series, when isolated in a very pure state, 
may be identified by the melting points, as is well known. The melting 
points for the series have been shown to lie on a smooth curve. In the 
case of mixed m-paraffins, Piper and co-workers (1931) showed that the 
purity of a paraffin sample, or the identities of the two or more paraffins 
in a mixture, cannot be determined from the melting point alone. 
However, the same investigators showed from experiments with -par- 
affins of known purity, and with known mixtures of these, that the purity 
of a sample and considerable information about the identities of com- 
ponents in an unknown mixture of n-paraffins can be determined by 
observing the difference between the transition temperature (on heating) 
and the melting point. It was also established by these same authors 
that there is an alternation in transition points (on heating), the differ- 
ences between the corresponding transition and melting points being 
from 1°C. to 2°C. greater for 7-paraffins with odd-numbered chains than 
for even-numbered. This alternation was observed through the u-paraf- 
fin series within the range from Cog to Cay. 

This alternation in the physical properties of the odd- and even-num- 
bered n-paraffins was cited by Miiller as confirmatory evidence support- 
ing his thesis of the zig-zag form of the molecule. From purely geometrical 
considerations, it would appear that a zig-zag molecular chain would be 
more in accord with the observed alternation in physical properties than 
would a spiral chain. It seems possible, however, that this alternation of 
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properties may be explainable on a chemical basis rather than geometri- 
cal. 

It is generally agreed that the transition point (on heating) marks the 
rupture of the inter-chain terminal bonds between adjacent end-to-end 
molecules. According to Miiller, this means that the chains retain their 
fixed relative positions in the crystal lattice, but are freed to rotate 
around their axes. It seems more likely that on reaching the transition 
temperature, the chains are not only freed to rotate on their axes, but 
also freed to move laterally in any direction normal to the c-axis, and thus 
to rearrange themselves from an orthorhombic packing into a hexagonal 
packing. When the temperature is raised to, and beyond the transition 
point, it is probable that the crystal would change from an orthorhombic 
to a hexagonal form, and the original @ angles of the crystal would 
change abruptly, decreasing toward the 60° angle as the melting point is 
approached. Observations supporting this interpretation will be described 
later. 

The n-paraffin molecule has generally been regarded as built up of a 
series of CH: groups, capped at each end by a CH; group. On this basis, 
there should be no alternation of physical properties of odd- and even- 
numbered members, except for the geometrical difference pointed out 
by Miiller. In a zig-zag chain which is even-numbered the C—C bonds at 
the ends of the chain are parallel; in the odd-numbered chain they are 
not. 

However, the writer has noted other apparently constant differences 
between odd- and even-numbered n-paraffins, which are not as readily 
explainable by the geometrical difference between the odd- and even- 
numbered molecule. Crystals of even-numbered n-paraffins have larger 
average @ angles than do crystals made up wholly or largely of odd-num- 
bered members of the -paraffin series. The even-numbered crystals also 
tend to be more strongly orthorhombic in their symmetry, and to have 
larger axial angles. These observations are only of qualitative value, and 
will be checked quantitatively when a sufficient number of pure ”-paraf- 
fins can be obtained for study. Nevertheless, differences in form, struc- 
ture, and crystallization habits do exist between odd- and even-numbered 
crystals, and these differences are more easily visualized as differences in 
the lattice dimensions than as differences of geometrical form of the 
molecules. 

It is suggested that the alternation of transition points on heating, as 
well as the above described differences in crystal form, which serve to dis- 
tinguish the odd from the even members of the #-paraffin series, are 
caused by chemical differences in the two types of molecule. If we regard 
the “building blocks” of the n-paraffin chain as C,H, and C.Hs groups, 
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rather than as the conventional CH, groups, then we have the basis for 
assuming a difference in chemical properties between the odd and even 
molecule. 

The evidence supporting this suggestion is indirect. Chibnall and Piper 
(1934) have shown rather conclusively that all long chain n-fatty acids 
and primary alcohols found in nature (various plant and insect waxes) 
have even-numbered chains, whereas all n-paraffin waxes from these 
natural sources have odd-numbered chains. These authors showed how 
the metabolic processes within plants might decarboxylate a fatty acid 
to form a corresponding n-paraffin with one less carbon atom. 

It is difficult to explain the synthesis by plants of even-numbered 
acids and alcohols exclusively unless the plants are assumed to utilize 
as “building blocks” short straight chain molecules of even number car- 
bon atoms. If we assume that in the plant metabolism long chain odd- 
numbered n-paraffin are synthesized, and some of them later carboxy- 
lated to form even-numbered acids, we are faced with the difficulty of 
explaining why both odd- and even-numbered paraffins are not found 
among the plant waxes. 

The assumption of primitive “‘building blocks’’ with even numbers of 
carbon atoms, perhaps synthesized during early metabolic processes in 
plants, would appear to explain the observations of Chibnall and Piper. 
It is suggested here that these primitive ‘“‘building blocks” probably in- 
cluded both C2H, and C4Hg groups. If only the latter group existed, then 
it would be difficult to explain why the fatty acids and primary alcohols 
with 26, 30, and 34 carbon atoms are about as abundant among the plant 
waxes as are those with 24, 28, and 32 carbon atoms. On the other hand, 
the C4Hgs group is the shortest chain which can define a right or a left 
handed spiral. It has been shown in the present paper that all the n- 
paraffins are dextrorotatory. If we assume that the primitive ‘building 
blocks”’ were C.H, groups exclusively, then it is necessary also to assume 
an inherent tendency in these groups to join in a left handed spiral. 

If the long chain molecules were all synthesized from original short 
chain groups, such as C;Hy, having an even number of carbon atoms, 
then an odd-numbered v-paraffin may be regarded as having at one end 
an incomplete or ruptured ‘‘building block” (such as a CH; group), and 
at the other end a normal C,H; (or CsHg) group. In a crystal lattice of 
an even-numbered u-paraffin, the inter-molecular bonds of end-to-end 
molecules would all be between like groups. In a lattice of odd-numbered 
chains, the end bonds would either all be between unlike groups, or be- 
tween CH; groups in one {001} interface, and between C,H; groups in 
the next interface, thus alternating throughout the crystal lattice. The 
latter arrangement would be in agreement with Miiller’s x-ray findings 
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that the c-spacing in odd-numbered »-paraffins is equivalent to the 
length of two molecules, and only equivalent to the length of one mole- 
cule in the even-numbered paraffins. The arrangement would also account 
for the observed alternation of transition points on heating, provided 
we assume that the CH;-CH; bond is less stable than the C,H;—-C2Hs 
bond. 


CHANGES IN CRYSTAL FORM ON HEATING 


Miiller (1930, 1932) measured the lattice dimensions of fifteen n- 
paraffins ranging from Cis to Cy at different temperatures by means of 
a series of x-ray photographs. He showed that the dimensions of the (a) 
and (6) axes of the unit cell changed with change of temperature. His 
results may be summarized as follows: 

1. On heating, the series Cz: to Cop changed from lower symmetry to hexagonal sym- 

metry, the latter being reached just below the melting point. 


2. On heating, the series Cig to C29, and those above C29 tended to approach hexagonal 
symmetry, but melted before arriving at hexagonal symmetry. 


3. Cx and C23 showed continuous change in (a) and (6) unit cell dimensions with in- 
creasing temperature up to the melting point. The m-paraffins Cos, C25, Cos, C27, Coo, 
C30, Cs1, Cs4, and Cy, showed abrupt transitions between room temperatures and their 
respective melting points. 


The unit cell dimensions of the »-paraffins at different temperatures 
as determined by Miiller may be used to compute the corresponding in- 
terfacial angles ¢ of the crystals. Some of Miiller’s results, thus converted 
to ¢ angles, are shown on the chart, Fig. 8. Also on the same chart is a 
curve for the ¢ angles of pure Octacosane (Cos) for different tempera- 
tures. This curve is based on direct measurements of the crystal on the 
microscope heating stage, and represents the mean of two series of ob- 
servations on two different microslides. The two original observed curves, 
however, are practically identical with the mean shown in Fig. 8. The 
errors of the observed points on the original curves were estimated to be 
less than +0.1° for the ¢ measurements, and approximately +0.2°C. 
for corresponding temperatures. 

In the first experiment, the crystal temperature was raised from 
24.5°C. to 58.5°C. over a period of 3.5 hours at the nearly constant rate of 
0.16°C. per minute. In the second experiment, the crystal temperature 
was raised from 23.5°C. to 56.9°C. over 2.25 hours at the average rate of 
0.25°C. per minute. 

Both original curves showed two definite transition points with rising 
temperature. The first (or lower) transition point was at 39.3°C. in each 
experiment. The second (or upper) and principal transition point was at 
57.3°C, in the first experiment, and at 56.6°C. in the second experiment. 
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The mean value, 56.9°C., is shown in Fig. 8. It will be observed that up 
to the temperature of 39.3°C., no change whatever took place in the 
crystal angles. From 39.3°C. to about 56.9°C. there was a linear decrease 
of about 1° in the ¢ angles. At 56.9°C. there was a very sudden decrease 
of about 3° in the ¢ angles. It will be noted in Fig. 8 that this ¢-tempera- 
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ture curve for Cos, based on direct measurements of the changing crystal 
angles, shows a close parallelism to the curve for C2:, computed from 
Miiller’s x-ray measurements of the unit cell. 

In the two experiments on Cos crystals, the temperature increase was 
carried to about 58°C. (3.5°C. below the melting point) and then al- 
lowed to decrease at uniform rates to room temperature over periods of 
two hours and four hours respectively. In each experiment the crystal 
{110} faces became somewhat irregular or warped at temperatures 
immediately above 56.9°C. (the upper transition, heating). During cool- 
ing below that point, still other irregularities appeared in the crystals, 
all of which had been completed on cooling to about 50°C. There were 
no further changes apparent in the crystals on cooling from 50°C. to 
room temperature. 

The changes during cooling were observed between crossed nicols, 
with gypsum plate. It was evident that the setting of the crystals was 
accompanied by an irregular reorientation of the lattice structure, and 
that the changes described were not merely surface or resorption phe- 
nomena. After cooling to room temperature, the measured crystals, as 
well as others on the same microslides, were studied conoscopically. The 
interference figures were for the most part diffuse, due evidently to the 
irregularities in the lattice structure. The most significant observation 
noted, however, was that the dextrorotatory optical activity was still 
present, though in many instances, the positions of extinction and of 
maximum double refraction had been shifted varying amounts with 
respect to the original axes of symmetry of the crystals. 


INTERPRETATION OF HEATING STAGE OBSERVATIONS 


The interpretations of the foregoing observations are summarized as 
follows: 


1. The first transition on heating (39.3°C. in Cog and about 40°C. in Cog, from Miil- 
ler’s data) marks a sudden “‘weakening”’ of the intermolecular end-bonds, permitting 
a progressive distortion of the orthorhombic lattice in the direction of hexagonal 
packing. 

2. The second transition (This corresponds to Piper’s “First Transition’”’—see Piper, 
1931) represents sudden rupture of the end-bonds of most of the molecules, and the 
initiation of a “mass movement”’ of the molecules toward hexagonal packing, as 
the ¢ angles of the crystal decrease rapidly toward the limiting value of 60°. If it 
were possible to bring all portions of the lattice to the transition temperature at ex- 
actly the same instant, the ideal curve would presumably be “flat’’ at this tempera- 
ture. In other words, the ¢ angles would decrease instantly to 60° when the crystal 
reached the transition temperature. It is also possible that the 60° ¢ angle could be 
reached by maintaining the crystal at exactly the transition temperature for a suffi- 
cient time to allow the molecules to move into the hexagonal pattern. 
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3. By allowing the crystal to cool before it had reached the melting point and before 
the entire lattice had reached the hexagonal form, a distorted lattice resulted, con- 
sisting of various segments, some orthorhombic, others more nearly hexagonal. This 
heterogeneity of the lattice is responsible for the observed warped and irregular crys- 
tal faces. 

4. The fact that the dextrorotatory behavior of the crystals persisted after the crystals 
had cooled from temperatures above the upper transition point, and in spite of the 
disarrangements in the lattice structure, confirms the interpretation that the optical 
activity of 2-paraffin crystals is produced by the internal structure of the molecule, 
and not by any peculiar orientation or arrangement of the molecules in the lattice. 

5. It will be noticed that Miiller’s curve for C2; (shown in Fig. 8) differs radically from 
the curves for C24, Cos, and C29. As mentioned above, Miiller found that C2; also 
behaved in the same manner as C23. Incomplete heating stage observations of mix 
n-paraffin crystals indicate behavior somewhat similar to Cz: and C23, which suggests 
that the samples of C2: and C2; used by Miiller may not have been pure compounds. 


GROWTH OF N-PARAFFIN CRYSTALS FROM SOLUTION 


One of the most interesting aspects of the study of normal paraffin 
crystals, is to be found in observing through the microscope their 
growth, either from solution at ordinary temperatures, or from the melt. 
The present paper, however, is concerned only with the crystals formed 
from solvent solution at temperatures well below melting and transition 
points. Under such conditions, many observations were carried out with 
different optical combinations, including plane polarized light, crossed 
nicols with parallel light (with and without gypsum plate), and cono- 
scopic observations of the interference figure. 

A complete summary of these observations would be unnecessarily 
lengthy, and for present purposes it is sufficient to record only the fol- 
lowing observed features, which are believed to have a bearing on the 
form of the crystal lattice, and the manner in which it is built up from a 
solvent solution of the paraffin. 


1. In oblique or reduced ordinary or plane polarized light, the crystals make their first 
appearance out of solution as full sized crystals; that is, the rhombic plate has at- 
tained its maximum dimensions in the plane of the a and b axes. The first appearance 
of the crystal is sudden, because at this stage of growth the crystal is a very thin 
rhombic plate, and is at the lower limit of visibility. The “‘suddenness”’ of appearance 
is therefore an optical illusion. During the course of minutes, the outline of the crys- 
tal becomes progressively clearer. The refractive indices of the crystal are in almost 
all cases higher than that of the surrounding liquid or amorphous medium. 

2. The same phenomena may be observed between crossed nicols, with and without 
a gypsum plate. With the gypsum plate, the interference colors, ranging between 
blue and yellow, are seen to increase gradually in intensities, the colors changing with 
the increase in thickness of each crystal. 

3. Conoscopically, the changes observed are similar. The interference figure, during 
early stages of crystal growth, is invisible, then very faint, and gradually strengthen- 
ing to normal visibility. During this change, the axial angle remains constant. 
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4. In the case of twinned crystals, the twins make their appearance simultaneously 
and in the same manner as described for the single crystal. 


The most probable interpretation of the phenomena described above 
is that the initial stage of the crystal growth is the attachment of mole- 
cules to the glass surface (both microslide and cover glass) as a single 
layer of molecules. The molecular chain axes are normal to the glass 
surface. The molecules arrange themselves in the orthorhombic packing 
found in the finished lattice. This primitive layer grows laterally until it 
attains the full (a) and (b) dimensions of the completed crystal. 

The second and subsequent stages of crystal growth consist in the 
consecutive additions of identical layers on the primitive layer, the 
process continuing until stopped by (a) interference of opposite glass 
surface, (b) interference from adjacent crystals, (c) exhaustion of paraf- 
fin in solution, or (d) by evaporation of the solvent and excessive viscos- 
ity of the surrounding solution. 


CRYSTALLIZATION UNDER OTHER CONDITIONS 


It is again emphasized that the foregoing description of the growth of 
n-paraffin crystals applies only to crystallizations from solvent solutions 
on glass microslides at normal room temperatures which are below the 
transition and melting points of the compound. The description also 
applies only to reasonably pure individual or mixed m-paraffins, free of 
oils or other organic materials which inhibit or modify normal growth of 
the paraffins. 

Under other conditions, the method of growth and the final crystal 
forms of the m-paraffins may be different from those we have described. 
For example, if during growth of a crystal the c-axis of the lattice is not 
parallel to the optical axis of the microscope, change in crystal dimen- 
sions may be visible during growth. When crystals are formed out of 
solvent solution at temperatures between the lower transition and melt- 
ing points, the lattice may be warped, so that the c-axis lies at angles 
other than vertical. A common habit under these conditions is the for- 
mation of ‘‘needles.’”’ One type of paraffin needle is formed by a rolling 
up or curling of the {001} layers around the axis of the needle. During 
growth, these needles are seen to increase rapidly in length and more 
slowly in width as successive {001} layers are added to the elongate, ci- 
gar-shaped body. The long axes of such needles usually lie horizontally 
on the microslide. The {001} layers are wrapped around the needle in a 
manner which may be compared to a long stack of conical parer drinking 
cups, or to the wrapper leaves on a cigar. As viewed through the micro- 
scope, the c-axis of the distorted layers of the lattice approaches the 
horizontal position toward the lateral margins of the needle, and the 
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vertical position toward the median or central portion of the needle. 
Thus the birefringence seen toward the two lateral margins of the needle 
is much stronger (y-@) than that seen along the central axis of the 
needle (8-a). As observed between crossed nicols, such needles present 
a characteristic appearance of two bright zones along the margins of the 
needle, separated by a darker zone extending along the needle axis. 
Needles of this type are seen in m-paraffins which have crystallized from 
the melt, and sometimes in those crystallized from solvent solution at 
ordinary temperatures, but with oil, asphalt, or some other organic ma- 
terials present in the solution as impurities. 

The formation of needles of the type described, and many other 
variations of crystal form and habit have been described in various pub- 
lished papers on the subject (see Heck, 1936; Kolvoort, 1938; Gray, 
1943; Rhodes et al., 1927; Tanaka et al., 1928; Ferris et al., 1931). These 
numerous other crystal forms are rarely or never observed when reasona- 
bly pure individual or mixed -paraffins are crystallized under a cover 
glass on a microslide at ordinary temperatures, employing such solvents 
as carbon tetrachloride, benzene, toluene, and xylene. 
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ABSTRACT 


The notable lack of quantitative information regarding errors in many mineralogic 
computations has prompted the writers to construct a series of diagrams from which such 
information can be readily obtained. These include various determinative methods for 
density, refractive index, birefringence, and optic angle Assuming a skilled operator, 
maximum values of the error in single observations are assigned to each procedure, so that 
the diagrams represent “worst” values of error, rather than a “probable” error as deter- 
mined from a series of observations. Comparison of the various procedures is based on this 
‘‘maximum expected” error. 
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INTRODUCTION 


It is a truism that all physical measurements are approximations 
which approach, in varying degree, an unattainable absolute value. 
As a result it is essential that the degree of approximation, or amount 
of fractional deviation, be adequately known. In the field of mineralogy 
everyone knows about this, but few have consistently done anything 
about it. For various reasons the subject as a whole has been treated in 
the manner of a poor relation, with much overstatement and some un- 
derstatement of the fractional deviations obtained (Hey, 1933). The 
present paper is an attempt to improve this situation. 

Deviations (also called uncertainties or errors) are classified as sys- 
tematic (periodic) or random (accidental). Systematic errors are in- 
herent in every experimental investigation (e.g. temperature variations 
in density and refractometric work) and can be controlled. This con- 
trol is in some cases fine enough so that this type of error is negligible 
compared with the random error. Random deviations arise through 
reading of scales, adjustment of extinction positions, etc. For a single 
observation the uncertainty of the measurement depends in part on the 
skill of the observer, in part on the character of the apparatus. If re- 
peated observations are made, the random error may be expressed quan- 
titatively by calculating a “probable error.” As there is a 50% probabil- 
ity of a measurement differing from the true value by more than this 
amount, the “probable error” of an observation may give a false im- 
pression of accuracy. In the present discussion ‘‘maximum expected 
error” is considered only—i.e., it represents the probable worst value 
which a skilled operator would obtain from a single measurement with 
a given apparatus. If it is known in advance how large this maximum 
expected error is likely to be and if this value, added to the known sys- 
tematic error, gives a total uncertainty which can be tolerated for the 
work in hand, the labor of computing a standard “probable error’? may 
be avoided. This maximum error is of course not a precise quantity, but 
the numerical values used in constructing the diagrams given here are 
believed to be large enough so that the curves give fractional deviations 
which would be exceeded only in rare instances. 

The general method followed is (1) differentiation of the basic formula 
with respect to all pertinent variables, (2) assumption of maximum ex- 
pected values for the differentials, (3) solution of the differential equation 
to determine the maximum expected random error, (4) calculation of the 
systematic error, (5) addition of random and systematic errors, (6) con- 
struction of a diagram. 

It is assumed that the reader is familiar with the determinative meth- 
ods under scrutiny and has acquired enough manipulative skill to avoid 
random errors appreciably greater than those shown in the diagrams. 
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In cases where the observer believes that the uncertainty in any of his 
own measurements differs markedly from the values used here, an ad- 
justed diagram can usually be made with little labor by referring to the 
proper differential equation. 


DETERMINATION OF DENSITY 


Four commonly used methods of density determination are consid- 
ered and their relative accuracies compared. The suspension, hydrostatic, 
and pycnometer procedures give relative density, or “‘specific gravity,”! 
whereas the x-radiation method measures density directly, without use of 
a reference fluid such as water or toluene. 

SUSPENSION METHOD. Although this method does not involve any 
computation and is limited in its application, its potential accuracy is 
worthy of note, especially where Clerici’s solution is used (Jahns, 1939). 
The random error in density inherent in Jahns’ refractive index-density 
curve is given as .01. A second random error arises from reading the 
Abbé refractometer, and may be assumed to be .0005. Using the slope of 
this curve to obtain the ratio of density and index change (dD/dn), the 
uncertainty in density is found to be .00044, negligible in comparison 
with the first random error. The total random error of .0104 is converted 
to a fractional deviation dD/D, this value being plotted against D in 
Fig. 2 (inset, top left). It is assumed that the systematic error arising 
from temperature variations can be adequately controlled. 

Extension of the density range through use of a glass float (of known 
weight and density) attached to the mineral fragment can be achieved, 
but with relative loss in accuracy. This case is not considered here (see 
Bannister and Hey, 1937). 

HyprosTATIC METHOD. A number of methods of this type have been 
described, of which two are discussed here. 

Jolly Balance. The expression for density is 


ie se 
wi [Wwe 


where Wa=weight of substance in air. 
Ww=weight of substance in water. 
Dw=density of water at room temperature relative to Dw=1 
at 4°Cs 


1 Although widely used, “specific gravity”’ is terminologically incorrect as a synonym 
for relative density, and is physically incorrect as a synonym for density determined by 
the x-radiation method. “Gravity” implies the weight of a body, i.e., the earth’s attraction 
for it, which is not an intrinsic property, whereas by definition density is the mass of a 
body per unit volume, the mass being an intrinsic and invariable property. 
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The fractional deviation is given by 
dD WadWw+ WwdWa i dDw 
D  Wa(Wa- Ww) Dw 


? 


where dWw and dW, are the uncertainties in determination of Wa and 
Ww, respectively, and dDw is the uncertainty in density of water at any 
temperature. As dW, and dWw are obviously equal they may be re- 
placed by dW. Furthermore, the uncertainty in Dw is relatively negligi- 
ble so that, finally, 

dD dW(Wa + Ww) 

D  Wa(Wa — Ww) 


which represents the fractional deviation arising from the random errors 
of reading the scale. The obvious source of systematic error would be 
temperature variations which would affect the density of the water. 
However, it is found that, between 15° and 25°C. (room temperature 
range), this error may be neglected. 

Figure 1 has been prepared from data obtained with the current ap- 
paratus offered by the Central Scientific Company. Two springs are 
provided to extend the range of the balance. The uncertainty of any in- 
dividual reading of the scale should not exceed one-quarter of a division. 
In order to compare directly with errors found in the other weighing 
methods, the springs were calibrated in grams, giving dW=.005 gm. for 
the light spring, and dW =.025 gm. for the heavy spring.” This conversion 
is of course not ordinarily necessary. Per cent fractional deviation is plot- 
ted against density for each spring, using four values of Wa known to be 
within the range of the scale. The sensitivity of each spring is considered 
constant within the range selected. It is seen from Fig. 1 that least ac- 
curacy (greatest fractional deviation) is obtained where Wa is low and D 
is high, and, conversely, that maximum accuracy is to be expected from 
the combination of high Wa and low D. Furthermore, the variation in 
dD/D with D is much less for high Wa than for low Wa. These observa- 
tions hold also, in varying degree, for the other balance methods dis- 
cussed (Figs. 2 and 3). 

Berman Balance. As with the Jolly balance the fractional deviation is 

dD dW(Wa + Ww) dDr 


Dyn Wal Wome Woes 


where dDz and Dy refer to toluene instead of water. At any given tem- 


? Figure 1 will probably require some adjustment for Jolly balances of different make. 
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perature there is an uncertainty dDr=.0005 in the density of toluene 
(Timmermans, 1912). The random error dW of a single scale reading on 


6 
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D 
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Fic. 1. Jolly balance. Variation in accuracy for different densities. Random error com- 
puted for Wa values of 1, 2, 3, 4 gm. (using light spring) and 10, 20, 30, 40 gm. (using heavy 
spring). Based on an error of .25 division in reading the scale. 


Fic. 2. Berman balance and (inset, top left) Clerici solution method (Jahns). Random 
error for Berman balance computed for Wa values of 5, 15, 20, 25, 50 mg. Based on dW 
=.01 mg. Full lines (15, 20, 25 mg.) cover the range of maximum sensitivity of the balance 
as specified by the manufacturer. Random error for Clerici solution method based on ex- 
perimental data given by Jahns (1939). 
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the balance should not exceed .01 mg. The systematic error arising from 
determinations made at various temperatures is relatively negligible in 
the range 15°-25°C. Surface tension deviation may also be neglected 
where the suspension wire diameter is a small fraction of millimeter. 
Figure 2 is therefore constructed directly from the equation given above. 
Optimum conditions of accuracy are realized only for Wa values from 
15 to 25 mg., as specified by the manufacturer. 

Where a determination is made with several small fragments instead 
of with a larger, single piece, accuracy will probably be lowered because 
of less perfect wetting. Unlike the pycnometric method, this systematic 
error can not be eliminated. 

PYCNOMETER METHOD. Density by the pycnometer procedure is given 
by 

(We — Wi)Dr 


Dp = 
We — Ws + DrV 


where W:= weight of pycnometer 
W.=weight of pycnometer+ powder 
W;= weight of pycnometer+powder-+ toluene 
Dr=density of toluene 
V=volume of pycnometer. 


If a silica-glass pycnometer is used, the uncertainty in V may be reduced 
by careful calibration to a relatively negligible value (Ellsworth, 1928). 
In the following expression for fractional deviation V is therefore consid- 
ered constant 


dD r, 2d W i —( W. — Ws; ) 
D We—-W, Dr \We — Ws + DrV/ 
The three uncertainties dW, dW2, dW; used in the first step of the dif- 


ferentiation are equal and have been replaced by dW in the final con- 
densed expression above.’ As the uncertainty in weighing on the ordinary 


5 As the writers have not succeeded in duplicating the expressions derived by one or 
two other authors for the pycnometric method, the derivation used here is given in detail. 


= (W2 — Wi)Dr ‘ 
W: — W; + DrV 
Taking logarithms of each side, 
log D = log Dy + log (Wz — Wi) — log (W2 — W; + DrV). 
Differentiating, 
dD fc dDr dW. — dWi tu dW. — dW; + VdDr 
D Dr’ W-W, W2— Ws+DrV 
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chemical balance should not exceed .0001 gm., this value is assigned to 
dW in constructing Fig. 3. dDr=.0005 as in Fig. 2. Room temperature 
variations produce relatively negligible systematic errors. The 10 cc. 
bottle described by Ellsworth (1928) is used to show variations in ac- 
curacy for powder samples weighing 1, 2, 5, and 10 gm.; the 0.1 cc. bottle 
described by Winchell (1938) is used similarly for powder weights of 0.1, 
0.2, and 0.5 gm. 

Every precaution must be taken to wet the powder thoroughly, as 
failure to remove all the imprisoned air will lead to deviations much 
greater than the random error from weighing. Ksanda and Merwin (1939) 
have described an excellent technique of this kind. 

Figure 3 shows that bottles of large capacity give smaller fractional 
deviations over the usual range of powder weights than do bottles of 
small volume, assuming equivalent conditions of weighing. This becomes 
obvious from inspection of the differential equation. Increase in V re- 
duces the dD7/Dr term, and increase in W.-W, (weight of powder in air) 
reduces the dW term. Offsetting this apparent superiority of the larger 
bottle, however, is the fact that its greater bulk of powder increases the 


Collecting the dW,, dW2, dWs, and dDy terms separately, 
dD Vv 
— = dDr{ — —- oar) - aw(———_ ) 
D Dr W.2 — W3 + DrV W2 — Wi 
1 
gapfial worsens ci caleoet be) 
Tie Wace pW t= We DV 


1 
+ eee Ww; + =) 


For determination of maximum error all the bracketed expressions must have the same sign. 
The dW, term is therefore made positive. (The minus signs within the dD and dW. 
brackets must be left unchanged, as Dy and W? occur in both numerator and denominator 
of the original density equation. Errors in their measurement would therefore be expected 
to compensate each other partially, as shown by the minus signs in the differentia] equa- 
tion.) Furthermore, as dW;:=dW2=dWs, we replace each by dW. Thus 


dD 


Bol pant 


Di(Ws — We + Drv) 
4 (W2 — Ws + DrV) + (W2 — Ws + DrV) — (We — Wi) + (W2 —W)) 
+ wl (W. — W)(W2 — Ws + DrV) 


=r ang Weep Lagoa 


DrdV 
If V is not considered constant the term [ W.LowseeiDrV 


| must be added totheabove. 
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difficulty of removing the last traces of trapped air. Homogeneity in com- 
position of the sample is also more difficult to attain with a large sample. 
It is probable therefore that in practice the accuracy of the two bottles 
will be about the same. 

In the micropycnometric procedure described by Bannister and Hey 
(1937) only 15 to 25 mg. of sample are required, with a reported accuracy 
of 0.5%. As part of this error stems from an uncertainty in the volume 
of liquid used, a greater error than that of Fig. 3 would be expected. 


Fic. 3. Pycnometer method. Variation in accuracy for different densities. Random error 
for 10 cc. bottle (Ellsworth type) and .1 cc. bottle (Winchell type). Full lines refer to the 
10 cc. bottle for mineral weights (in air) of 1, 2, 5, 10 gm.; broken lines refer to the .1 cc. 
bottle for mineral weights (in air) of .1, .2, .5 gm. Based on dW=.0001 gm. and dDr= 
.0005. 


However, uncertainties regarding the homogeneity and wetting of the 
powder may be better controlled in the Bannister-Hey method than in 
the others, and too much emphasis should not be placed therefore on the 
purely mathematical aspects of the case. 

X-RADIATION METHOD. Density by this method is expressed by 


MZ 
Dise-as 
AV 
where M=molecular weight of the material 
Z=number of formula-units in a unit cell 
V=volume of unit cell 


A= Avogadro’s number (6.02 X 1078). 
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As V is the single variable, 
AD AV 


D V 


In a cubic crystal, a9 =d(h?+h2+1?) 2 
where a)= length of side of unit cell, 
d=interplanar spacing of planes with indices A&l. 


Aa Ad A 3 
Therefore —- = —— But aN = as) = A) = ee 
a0 d V ao? ao d 


From Bragg’s equation d=)/2 sin, 
where \= wave length of x-radiation used, 
6= diffraction angle. 


The fractional deviation of d is Ad/d=—cot 6A. Expressing the frac- 
tional deviation in density in terms of the diffraction angle, AD/D =3 
cot 646 (for cubic crystals only). 

Application of this equation is confined here to powder photography. 
The error in reading a millimeter scale placed on a line of the film should 
not exceed .25 mm. for any single measurement. Therefore the uncer- 
tainty in the distance between such a line and its mate on the other side 
of the line of zero diffraction is 0.5 mm. As this total distance in a camera 
of 57.3 mm. diameter is equivalent to 46, the error in @ is .125°. That is, 
A@= .125/57.3 =.00128 radians. The curve R for random error in Fig. 4 
is based on this value. 

The four systematic errors inherent in x-ray crystal diffraction work 
—absorption of radiation by the specimen, eccentricity of the sample with 
respect to the camera, film shrinkage after development, and uncertainty 
in the camera radius—must each be considered on its merits. As details 
are given by Buerger (1942) for such corrections the matter is taken up 
only briefly here. Use of a two-hole film (Straumanis method) permits the 
camera radius and film shrinkage errors to be combined, and for a care- 
fully constructed camera and film dried at room temperature, the error 
from this source is relatively negligible. It does not appear in Fig. 4. 

The eccentricity error is small but appreciable. Its curve (E in Fig. 4) 
is based on a maximum expected displacement of the sample relative to 
the true center of the camera of .025 mm. The equation used is given by 
Buerger (1942, p. 413). The deviation may be positive or negative. 


4 To avoid confusion with the interplanar spacing symbol d, A replaces here the custom- 
ary differential symbol. 
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The absorption error (Buerger, 1942, p. 414) is based on complete 
absorption, a specimen diameter of .25 mm., a pin-hole—specimen dis- 
tance of 10.6 mm., and camera diameter 57.3 mm. Its curve (A in Fig. 
4) shows an order of magnitude for AD/D similar to the curve R, and of 
the same sign. 

The total value of the maximum expected error is shown by curve T 
and has been drawn using positive values for the individual errors. As 
the random error may be opposite in sign to the absorption error, and is 


Fic.4. X-ray method. Variation in accuracy for different diffraction angles. Random error 
(curve R) based on an uncertainty of .25 mm. in reading the position of a diffraction line. 
Absorption error (curve A) from data in Buerger (1942). Eccentricity error (curve E) like- 
wise from data in Buerger (1942). Curve T shows total error where R, A, E are assumed to 
have the same sign. 


of the same order of magnitude, the total error may, however, be very 
small. Curve T represents the worst possible case. The diagram does 
not include values smaller than 50°, as the accuracy decreases very rap- 
idly in this range. 

COMPARATIVE ACCURACY. Excluding the x-ray method, comparison of 
the remaining procedures shows that the pycnometer method has the 
greatest inherent accuracy, the Jolly balance method the least. Choice 
of method is dictated by the nature of the material and the degree of ac- 
curacy required. For most work the accuracy range of the Berman bal- 
ance is adequate, provided that a single fragment of suitable weight can 
be obtained. 

Direct comparison of these weighing methods with the x-radiation 
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method is not possible, as the latter gives the density of the unit cell 
whereas the former give values for aggregates of unit cells. The individual 
cells in these aggregates are rarely, if ever, packed together perfectly, 
with the result that lower density values are obtained than from the x-ray 
procedure. If it is desired to contrast the two density values for a given 
material, a weighing method should be selected which gives an accuracy 
comparable with that indicated by the @ value chosen in Fig. 4. In gen- 
eral, the Berman balance would satisfy this requirement. 


DETERMINATION OF REFRACTIVE INDEX 


As this paper is concerned with errors in computation only, the ac- 
curacy of direct-reading instruments will not be discussed. In the case of 
refractometers the accuracy of the widely used Abbé type has recently 
been reviewed by Tilton (1942) in an informative paper. It need only be 
stated here that the ordinary Abbé used in mineralogical laboratories 
reads correctly to +.0002 in its lower range, but is considerably less 
accurate in its upper range. Other direct-reading refractometers (Fisher, 
Nichols, etc.) are usually reliable to +.001. 

MINIMUM DEVIATION METHOD. In the field of precise refractometry the 
minimum deviation method is preferable to any other. Where moderate 
precision is required, as in standardization of oils or melts for use in the 
immersion method of refractive index determination, it affords the only 
universal method for the entire range of indices. The accuracy of the 
method has been discussed by Tilton (1935) who lists 19 sources of error 
for determinations of the highest precision (+ .000001). As mineralogists 
seldom require this accuracy for crystal measurement, and never for im- 
mersion oils or melts, only the random error of reading prism and mini- 
mum deviation angles need be considered here. 


The index of refraction is given by 


Ayan 


where A=angle of the prism. 
Dm= Minimum deviation angle. 


Differentiating partially with respect to A, dn=—1/2 sin Dm/2 cosec’? 
A/2 dA, and partially with respect to Dm, dn=1/2 (cos Dn/2 cot A/2 
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—sin Dn/2) dDm. As it is probable’ that dA=dDm, we replace each by 
dr, so that the total random error is 


1 A A+ Dn oe ne 
Z(dn) = =| eosee — cos ————— — cosec? — sin —— | dr. 
2 2 Ve 2 2 


In constructing Fig. 5 the error in the prism and minimum deviation 
angles is assumed to be +1’, so that dr=.000291 radians. The second 
term in the differential equation is considered positive in order to give 


0-0016 


0-0000 
1:3 17 21 ‘ 25 29 
Fic. 5. Minimum deviation method. Variation of random error in for a series of n values. 
Computations shown for prism angles of 25°, 30°, 40°, 50°, 60°, with the % accuracy given 
for each. Steep curve transecting the prism curves is the locus of limiting values of prism 
angles for determination of a given refractive index. Diagram based on an uncertainty of 
+1’ in determination of prism and minimum deviation angles. 


dn its maximum expected value.® Most single determinations will have 
greater accuracy than that indicated by the graph. In working with crys- 
tal prisms the deviation is likely to be much less than 1’, but with the 
ordinary hollow cell used for oil determinations this is unlikely, and for 
crude prisms composed of S-Se mixtures may be much greater than 1’. 
Figure 5 shows the total maximum expected error for prisms of varying 


5 This may not be true if the prism angle is determined from readings of the table circle, 
and the minimum deviation angle from readings of the telescore circle. The table angles of 
small spectrometers and goniometers tend to be less accurate than telescope angles. How- 
ever, if the prism angle is determined by the split-beam method, no table angle need be 
used, and except in high-precision work gives very satisfactory results (Guild, 1923). 

® The small numerical range of the common refractive indices makes it unnecessary to 
plot the variation in fractional deviation for different values of n. 
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angle over a wide range of indices. For a given index greatest accuracy 
(least dn/n%) is attained where large prisms are used. Inspection of the 


formula 
i nee fs A 
n = sin ————— / sin — 
2 2 


indicates a maximum value of n for any given prism angle. The boundary 
curve in Fig. 5 shows this graphically. For example, the limiting refractive 
index determinable with a 60° prism angle is 2.0. In practice, since the 


00005 


Fic. 6. Minimum deviation method where prism angle is 60°. Curve P shows random error 
arising from prism angle uncertainty of 1’. Curve Dy shows random error arising from 
minimum deviation uncertainty of 1’. Upper curve shows the total error where P and Dy 
are assumed to have the same sign. 


proportion of transmitted light at grazing incidence is greatly reduced 
relative to the reflected portion, it is never possible to proceed to the limit 
shown in the diagram. 

Figure 6 gives greater detail for a 60° prism. Curve P shows the 
variation in prism angle error for various indices; Dm shows dn for dif- 
ferent indices arising from measurement of the minimum deviation an- 
gle. At the limiting index of 2.0 the prism angle error reaches itsmaximum 
whereas the minimum deviation error becomes zero. The upper, straight 
line represents the total maximum expected error. If the Gifford method’ 

7 This method (Gifford, 1902) makes use of a closed prism whose angles approximate 
60°. Minimum deviation angles are measured across each prism and an average computed. 
Angle A may then be assumed equal to 60° without any accurate measurement of the in- 
dividual prisms. An uncertainty in A of the order of magnitude of 1° causes an uncertainty 


of only +.0001 in m. As there need be no difficulty in maintaining the prism deviations at 
much less than 1°, the resultant prism angle error using this method becomes negligible. 
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is used, the error in prism angle becomes negligible and the total uncer- 
tainty is shown by Dm. In all cases where the demands of this method 
can be met, increased accuracy will result, particularly for refractive in- 
dices approaching 2.0 in magnitude. 

Where the goal of accuracy is in the neighborhood of + .0001, system- 
atic errors become significant only for determinations made on liquids. 
Temperature corrections are necessary, particularly for methylene iodide 


2 
13 14 is, 16 17 18 


Fic. 7. Grazing incidence method (Pulfrich refractometer). Variation in dn with n shown 
for two prisms of refractive index as marked. Approximate % accuracy at curve maxima 
is also given. Based on an uncertainty of +1’ in reading the scale. 


and mixtures containing it, and require control to at least 0.1°C. if the in- 
dex error limit of .0001 is to be maintained. In addition the hollow cell 
must be constructed of optically flat glass. Suitable material may be se- 
lected according to the method suggested by Larsen and Berman (1934). 
In Figs. 5 and 6 it is assumed that these systematic errors can be reduced 
to negligible proportions and that only random errors remain. 

GRAZING INCIDENCE METHOD. This method is incorporated in the 
widely used Pulfrich refractometer. Although less used by mineralogists 
than by chemists (Gibb, 1942, p. 333), and inherently less accurate than 
instruments used for precise minimum deviation work, a brief compari- 
son is not out of place. Refractive index is given by 


n= (np? — sin? r)¥/? 


where n,=refractive index of prism, 
r=angle read on divided circle. 
Differentiation gives 


| sin 2r 
dn = — | ————————— Idr 
2(np? — sin? a 
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The error in a single reading in the middle range of r values is not likely 
to exceed 1’, so that dr=.000291 radians. Figure 7 shows a graphical 
solution of the differential equation for two prisms of n=1.620 and n 
= 1.752. For each prism there is an r value which gives maximum frac- 
tional deviation (.07% and .06% as indicated). Thus for material of 
n=1.6 the maximum expected error would be only about one-half as 
great where the lower index prism is used as where the higher one is used. 
For high and low r values, however, dr may easily exceed the value used 
in the diagram, so that the accuracy indicated in the lower part 
(dn <.00004) may not be realized. However, even the worst values at the 
noses of the curves are as good, or better, than the accuracy to be ex- 
pected from ordinary four-figure minimum deviation measurements. 

UNIVERSAL STAGE EXTRAPOLATION METHOD. In cases where one critical 
index alone can be determined in a given oriented optic symmetry plane, 
measurement of an intermediate index at a known angle in this plane 
makes possible the computation of the other critical index. This method, 
first used by Pauly (Johannsen, 1918, p. 266), and later developed by 
Emmons (1943, pl. 10) depends on the following equation 


Ne(n;? sin? @ + nx? cos? ¢)!/? 
tx ) 
Nh 


where nx= unknown critical index, 
n, = known critical index, 
n2= intermediate index, 
¢=rotation angle used in obtaining no. 


Assuming that the uncertainties in n; and nz are the same, and using 
dn to represent both, the differential equation is: 


xe 2 il J 
dns [— Ps ee # -)| ance (1 a =~) cot odd. 
Nx Ne No? Ne nh ny’ 


If dn=.001 and the uncertainty® in ¢ is 1° (dg=.0174 radians), the 
equation appears graphically as in Fig. 8. It is assumed that ni>nx. 
In Fig. 8 dnx is plotted against (ni-nx) for different values of $ using 
the computations made for n}=1.6. For other values of ni no appre- 
ciable difference in the deviations is obtained. For ¢=90° inspection 
of the differential equation shows that dnx=(nx/n2) dn and since ne 
now equals nx, then dn, =dn=.001 as shown in Fig. 8. The total dn term 
remains fairly constant at .001+, most of the variation being shown by 


8 The uncertainty in @ stems from the initial random error in orienting the optic sym- 
metry plane rather than from any inherent defect in the OEW axis of the stage. 
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the d¢ term. In general, dn, is least where large angles are combined with 
low partial birefringence. Except for angles of less than 20°, however, 


dn, changes little where partial birefringence is low. 


009 


‘007 
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Fic. 8. Extrapolation method. Variation in dnx for different partial birefringence (mi 


—nx). Shown for three values of ¢, where n1=1.6. Based on dn=.001 and an uncertainty 
in ¢ of 1°. 


Fic. 9. Optic angle method. Variation in dy (or da) for different values of y—6 (or a—B). 
Four selected values of the optic angle are used. Diagram based on 8=1.50 and an optically 
positive crystal. Full lines refer to da only; broken lines to dy. The uncertainty in deter- 
mination of the known refractive indices is .001 and in the optic angle 1°. 
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If ny<nx, dnx increases by a few per cent for high nx and birefringence 
values; for low values the change is inappreciable. No diagrams for this 
case have been constructed. 

OPTIC ANGLE METHOD. Mertie’s nomogram (Mertie, 1942) showing the 
relations between a, 8, y, and V, based on ellipsoid geometry, is a neces- 
sary prerequisite to the deviation diagrams of this paper. As differentia- 
tion of the full equation gives a long, involved expression, the approxi- 
mate equations 


B—a 
sin? V = ——— for optically positive substances, 
oY te Oe 
; 0 
sin? V = 
yY—a 


for optically negative substances, 


are used here. Differentiation of these with respect to a, 6, y in turn, 
gives 


@ (8 — vy) sin 2VdV 


da = dn 
cost V F a = " 
or optica ositive substances. 
ty a ang Osim VV ie ea iti 
RA ee 
sint V 
— y) sin2VdV 
eae (6 ” 
sint V 5 5 
for optically negative substances. 
i sneer (8 — a) sin 2VdV 
— n —— ee 
fa cos! V 
dg = dn + (y — a) sin 2VdV for either positive or negative substances, 


where dn substitutes for d8 and dy, d@ and da, da and dy, respectively, 
and is assumed equal to .001. The uncertainty in direct measurement of 
2V need not exceed +1° if both optic axes are found. If only one axis can 
be oriented the uncertainty will be + 2°. For low to moderate refringence 
and birefringence 


da(+) =dy(—) and dy(+) = da(—). 


Figure 9 is based on a constant value of 8=1.50, an optically positive 
substance and dV=.0174 radians (uncertainty in V of +1°). It shows 
the variation in da with change in a-8, and the variation in dy with 
change in y-8, for selected values of V. As would be expected, the closer 
a (or y) approaches £ the smaller is da (or dy). The rate of increase in 
da (or dy) is greater for small V values than for large. The percent frac- 
tional deviation (accuracy) follows the same trend. Figure 10 is similar 
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to 9, but is based on B=1.90. The same type of variation persists, but 
da and dy are on the whole several times larger. For 6 values other than 
1.50 and 1.90 interpolation or extrapolation will give a reasonably good 
result. 


0:06 0-08 


Fic. 10. Optic angle method. Same as Figure 9, but based on B=1.90. 


Fic. 11. Optic angle method. Variation in d@ for different values of total birefringence 
(y—a). Four selected values of the optic angle are used. Based on uncertainties of .001 in 
determination of the known refractive indices, and of 1° in the optic angle. 


Figure 11 is of similar form, but the variation is simpler in character. 
The curves are valid for any value of 8, and for both positive and nega- 
tive substances. In contrast with Figs. 9 and 10, d@ is smaller for small V 
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values than for large. The overall range of df is also considerably less than 
for da and dy. 

If the refractive indices are determined by precision methods, dn will 
be much less than .001. However, as the immersion method is rapid and 
widely used, the diagrams have been constructed for dn=.001. 

COMPARATIVE ACCURACY. It has already been pointed out that the min- 
imum deviation method is capable of greater precision than any other. 
However, the small spectrometers used for this method in most labora- 
tories give results which are usually less accurate than those obtained 
with the Pulfrich refractometer, although adequate for practically all 
mineralogic work. 

In contrast, the extrapolation and optic angle methods, insofar as they 
depend on immersion in oils, are much less accurate, giving third decimal 
place accuracy only. 


DETERMINATION OF BIREFRINGENCE 


GENERAL COMPENSATOR METHOD (Berek). Determination of refractive 
indices to +.001 (immersion method) yields birefringence values accu- 
rate to +.002. This direct method may be compared with an indirect 
method which makes use of the Berek compensator. Birefringence is de- 
termined by: 

1p 


t/cos 6’ 


where I’ =retardation in the unknown for a particular orientation. 

t =thickness of thin section. 

6’ =angle between path of light and normal to the thin section. 
By differentiation, dB/B=dI’/I’—tané’d 6’—dt/t. The thickness 
t is found from a reference mineral of constant B (e.g. quartz), and 
t=I"’’ cos 6’’/B, where I’ and 6” are analogous to I” and 6’. Differentiat- 
ing, dt/t=dI’’/I’’—tan 6d 6’. By substitution in the original dB/B 
equation, and assuming positive values for each term (to give the maxi- 


mum deviation) 
dB dr’ ake 


B = Tr + tan 6'd 0’ + py + tan 6 'd.0% 


A general expression for dI’/T’ and dI’’/I'” may be obtained from the 
basic compensator formula, log '=log C+log f (i) where C=constant 
of the compensator;? f(i) =function of angle of rotation i of compensator 
drum. As f(i) =sin?i+.204 sin‘i+ --- for a calcite compensator (Berek, 
1913), differentiation of the compensator formula gives: 


9 This constant must be determined for each instrument. 
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dr sin 2i(1 + .408 sin? i)di 


ip: sin? i + .204 sin‘ i 


The maximum expected uncertainty of a single reading of the com- 
pensator drum is assumed to be 0.2°. The uncertainty in 2i would then be 
0.4° and for i would be 0.2°. Then di=0.2/57.3 =.00349 radians. Curve 
dI'/I in Fig. 12 shows the per cent variation in fractional deviation with 
change in I’.!° It is seen that accuracy is much increased where T is not too 
small. 


oy 
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Fic. 12. Berek compensator. Variation in dl’, and % accuracy of I’, for different values 
of I’. Based on an uncertainty of +.2° in reading the compensator drum. 


The uncertainties in 6’ and 6”’ may be assumed equivalent. Their mag- 
nitude depends on (1) error made in orienting the sections, (2) error in 
obcaining 6’ and 6” on the stereographic net, (3) error in the final hemi- 
sphere correction. Error 3 applies only to the unknown, as the refractive 
indices of the reference mineral are known. For the worst case, therefore, 
a maximum uncertainty of +2° is not improbable. Where 6’ and 6” are 
both small, however, errors 2 and 3 will be negligible and the uncer- 
tainty should not exceed + 1°. The general importance of small 6 values is 
indicated clearly in Fig. 13. 

No simple diagram will show the variation in accuracy of B, but the 
total fractional deviation in any given case can be found by adding the 
values of the four terms in the equation for dB/B. 


10 The values in Fig. 10 are about double those given by Berek (1924, Fig. 22, p.46). 
His calculation, however, is of a “mean” error, not a maximum expected deviation. 
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The preceding discussion implies strict uniformity in thin section 
thickness. As this condition is seldom realized, and, as it is physically 
impossible for the reference mineral and the unknown to occupy the 
same area of the section, a systematic error is thereby introduced which 
has not yet been considered. In the case of a strongly wedge-shaped sec- 
tion, detection by means of interference colors is easy, and the mount 
discarded if necessary. The thickness of apparently uniform sections can 
be rapidly checked by orienting three or four widely separated reference 
grains and noting the discrepancy in thickness obtained. Where this is 
small, e.g., .002 mm. across the whole section, it can be shown that, if 
unknown and reference mineral lie reasonably near to each other, a neg- 
ligible error in B will result. In the case of mineral fragment mounts a 
reference mineral can always be added, so that close spatial association 
with the unknown is realized without difficulty. 

In summary the accuracy of the indirect method increases where (1) 
thick rather than thin sections are used, (2) measurements are made on 
grains requiring the least angular rotation for their orientation, (3) the 
reference grain is as close as possible to the unknown. 

EMMONS EXTRAPOLATION METHOD. A new procedure described by Em- 
mons (1943, plate 13), based on the Biot-Neumann equation, permits 
determination of total birefringence in any crystal where the optic plane 
is vertical. The compensator used may be of fixed retardation (standard 
gypsum or mica plate) or variable (Berek, graduated wedge, etc.). As 
the analytical work on which Emmons’ plate 13 is based has not been 
published, no diagram of maximum error appears in this paper. Com- 
ment on the probable accuracy of the method will be found in the follow- 
ing section. 

Optic ANGLE-COMPENSATOR METHOD. The approximate relation 


Bp 
By= ———_——— 
sin 6 sin 6’ 
developed by Biot and Neumann (Johannsen 1918, p. 351), where Br 
=total birefringence of the crystal, Bp=partial birefringence (as meas- 
ured on a randomly oriented crystal section), @ and 6’ are the angles sub- 
tended by each optic axis with the normal to the crystal section, becomes 
By=(Bp/sin? V) where the optic plane is parallel to the microscope 
axis, with X or Z vertical, thus making 2=V=0’. A graphical solution 
of this equation is given by Emmons (1943, plate 11). Where Br>0.1 
there will be an appreciable correction to apply to V (see Larsen and 
Berman, 1934, Fig. 1, or Emmons, 1943, plate 13). By differentiation, 
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dBr dBp 
= — 2 cot VdV. 
Br Bp 


dBp/Bp consists of four terms as described for the general compensator 
method (Figs. 12 and 13). The additional term 2cot VdV is shown graph- 
ically in Fig. 13 for two values of dV. The uncertainty in measuring an 
optic axis position need not exceed + 1°; therefore if both axes can be 
measured directly, the error in 2V is + 2° and in V is +1°. If, however, 
only one axis position can be measured, the uncertainty in V is + 2°. Both 
these cases are included in Fig. 13. It is obvious that 2cotVdV is least 
where V is large and that the method fails badly for small V values. 


30 


2 cotVdV 


0° 10° 20° 30° 40° 


Fic. 13. Trigonometric terms used in the compensator methods. Tan 6d@ curve (general 
compensator method) based on an error of +2° in @. 2 cot VdV curves (optic angle-com- 
pensator method) based on an uncertainty in V of +1°in curve (a) and of +2°in curve (b). 


COMPARATIVE ACCURACY. Inspection of the differential equations de- 
rived for the first two birefringence procedures indicates that the general 
method is potentially superior to the optic angle method. Presence of 
the 2cotVdV term in the latter decreases the relative accuracy mark- 
edly. Total birefringence may be more accurately determined therefore 
on an optic normal section using the general method. 

It is probable that the fractional deviations in the Emmons extrapola- 
tion method are of the same order of magnitude as those in the optic an- 
gle method, i.e., somewhat less accuracy is to be expected than in the 
general method. 


Comparison of the general method (indirect) with direct determina- 
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tions of B may be made in an approximate way as follows. Where re- 
fractive indices are determined to +.001 (immersion method) the un- 
certainty in birefringence would be +.002. For B=.01 the fractional 
deviation is then 20%; for B=.10 it is 2%. Inspection of the range of 
values in Figs. 12 and 13 for the indirect general method indicates that 
a total fractional deviation in birefringence of 15% might commonly 
occur. By using thick sections and selecting grains which do not re- 
quire appreciable rotation for their orientation, the above value might 
easily be reduced to 5% as a maximum expected error. Comparison of 
the two methods depends therefore on auxiliary data, particularly on the 
magnitude of B and @. For low birefringence the indirect general method 
may, in many cases, be preferable; for high birefringence the direct 
method is superior. Values of B obtained from precise prism determina- 
tion of refractive indices would, of course, be superior in all cases to the 
indirect method. 


DETERMINATION OF Optic ANGLE 


Accurate information regarding the size of the optic angle in biaxial 
crystals is less important than accurate determination of refractive in- 
dices. However, as a number of procedures for optic angle determination 
are in common use, it is advisable to consider their relative merits. 

COMPENSATOR METHOD. The approximate relation between optic angle 
and birefringence, sin? V= Bp/Br is equivalent to sin ?, V=I'p/I'r, where 
retardations are computed for a common thickness of section, I'p is the 
retardation in a section normal to X or Z, and Ir is the retardation ina 
section normal to Y. Wright’s graphical solution of this equation is given 
by Emmons (1943, plate 11). Differentiating,” 


T'rdI'p = TpdI'r 
pk T’y sin 2V 


dV 


Using the Berek compensator, dI'p and dI'r are obtained from Fig. 12. 
In Fig. 14 dV/V (%) is plotted against I'v for several selected values of 
V. In general, accuracy is seen to be greatest for large retardations; in 
particular, where V is large as well. 

As the basic formula above is an approximation, an appreciable cor- 
rection to V must be applied where B>0.1 (Larsen and Berman, 1934, 
Fig. 1). This correction, however, would result in negligible change in 
Fig. 14. 

REFRACTIVE INDEX METHOD. The approximate equations cos? V4 
= (B-a)/(y-a) =sin? Vy, and sin? Va=(y-8)/(y-a) =cos’ Vy may be used 


11 The negative sign in the expression is changed to positive for computation of maxi- 


mum error. 
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to determine V if all three indices are known. The subscript a refers to a 
negative, to a positive mineral. Differentiating cos? Va, 


(y — a)dB + (6 — y)da + (@ — B)dy 
(y — a)? sin 2V 


dV = 


0 400 800 _ 1200 1600 2000 
ro 


Fic. 14. Compensator method. Variation in accuracy of V for different values of total 
retardation (I'y). Diagram shows curves for four values of V. 


Asa<B<vy, dV will have its maximum value where df is positive, da and 
dy are negative. The equation is then rewritten 


_ (yo a)dpe BE dee Baye 
= (y — a)? sin 2V 


dV 


Furthermore, da=d8=dy, permitting the substitution of dn, so that 


2dn 57.3 X 2dit 
or dV (degrees) = 


dV (radi = ——____-— Seared eames 
(radians) (pee a)? sin 2V (y - a)? sin 2V 


The uncertainty in n for the immersion method need not exceed 
+ .001, and for precise prism work +.0001. Figure 15 is constructed for 
these values of dn, the curves representing various values of V as 
marked. In general, accuracy is greatest for high birefringence; in par- 
ticular for high V values. It is obvious that an optic angle computation 
where dn=.001 is useful for relatively few substances (see Mertie, 1942), 
and that prism refractometry is required if any precision in V is expected 
by this method. Even with precision methods the accuracy is low if V is 
small. 
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EXTINCTION ANGLE METHOD. This method, first made practical by 
Berek (1924) and later amplified by Dodge (1934), provides a means of 
determining the optic angle independently of either refractive indices or 
retardation. It is discussed in two parts as outlined below. 

Optic plane horizontal. This is the most useful case. The geometrical 
relations are as follows (Dodge, 1934): 


K + K’ 
= ———— tan K = tan (6 + V) cosx, tan K’ = tan ($6 — V) cos x, 
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Fic. 15. Refractive index method. Variation in accuracy of V for different values of total 
birefringence (y—a). Curves shown for optic angle values of 10°, 20°, 30°, 45°. Based 
on refractive index uncertainties of .001 and .0001. 


where the following rotations are made, in the order listed— 

¢@ is a fixed angle of rotation on the outer vertical axis of the stage; 

x is a fixed angle of rotation on the outer east-west axis of the stage; 

¢ is an angle of rotation to extinction made on the microscope axis. 
The graphical solution of these equations will be found in Dodge (1934) 
and also in Emmons (1943). 

For the so-called ‘‘normal’’ case, where $=45°, and x is set equal to 
some given angle, we assume V and ¢ to be the only variables” in the 
above relation. Differentiating, 


2d¢ 


dea a eres 
sin 4¢ tan 2V 

12 The uncertainty in x arising from incorrect orientation of the optic symmetry plane 
is lessened very considerably here by the large ¢ value used, and is neglected in the differ- 
ential equation. Where ¢ is small x cannot be considered constant. 
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Fig. 16, showing dV/V plotted against V, contains two curves (full lines) 
derived from this equation, for values of x=54.7° and 70°. The uncer- 
tainty in ¢ need not exceed +1° if the optic axes are more nearly hori- 
zontal than vertical for the orientation in which ¢ is determined. If, how- 
ever, the opposite condition holds, as it will in many cases, the uncer- 
tainty in ¢ is likely to be more nearly +2°. As Fig. 16 is based on this 
larger value, dV/V may be halved if the optic axes approximate a hori- 
zontal position under the conditions just stated. It is seen that high V and 
high x increase the potential accuracy of the measurement. Two incom- 


30 
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Fic. 16. Extinction angle method (optic plane horizontal). Variation in accuracy of V for 
different values of V. Full lines for 6=45° and x=54.7 and 70° as marked. Broken lines 
for x=70° and @=15° and 30°. These four curves based on an uncertainty in ¢ of +2°. 
Curves D, and Dz, for direct readings of optic angle, are based on uncertainties in the optic 
angle of 1° and 2°, respectively. 


plete curves for ¢=15° and 30° (both for x=70°) have also been plotted, 
using the slopes of the proper curves given by Dodge (1934, plates I and 
IT). Although these curves indicate higher accuracy than where 6=45°, it 
must be recalled that x is not constant for such small values of ¢ (see 
footnote to previous paragraph). As the differential equation for 6<45° 
and x variable is extremely cumbersome, no computations have been 
made. However, the discrepancy in fractional deviation between curves 
based on $= 45° and $= 15° is less than Fig. 16 indicates. 

Figure 16 includes also two curves, marked D, and Dz, which show the 
variation in accuracy where V is measured directly. D, is based on the 
case where both optic axes are measured and the uncertainty may attain 
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a value of + 1°; Ds is for the case of one measured axis and an uncertainty 
of +2°. Assuming an uncertainty of 2° in ¢, and $= 45°, it is seen that 
the computed values of V approximate the accuracy of the direct method 
only for the highest values of V; for small optic angles a small value of J 
is required to attain comparable accuracy. 

Optic Plane Vertical. This orientation is considered here only for the 
case of a horizontal acute bisectrix. Where the acute bisectrix is vertical 


10 20. -y~~S(80 40 


Fic. 17. Extinction angle method (optic plane vertical). Similar to Fig. 16. 


it is usually possible to measure both axes directly and thus avoid indi- 
rect methods. The basic equation (Dodge, 1934) is: 


cos 2¢ cos x 


sin? x sin? x 


1/2 

sec Ve— (sin’ o+ sin 26 cot 26) 

where @, x, and ¢ are the same quantities used in the previous section and 
rotations are made in the same order. As before, ¢ and V are assumed 
to be the only variables, although this is true only where x is not appre- 
ciably affected by incorrect setting of the optic symmetry plane. Dif- 
ferentiating, 

cos x sin 26 cos Vd@ 


cos 2¢ cosx 


dV=— 
sin? x sin? 2¢ tan V (sin’ o+ 


1/2 
sin 26 cot 26) 


sin? x sin? x 


For selected values of ¢ and x this expression becomes much less formida- 
ble in appearance. The uncertainty in ¢ is taken as + 2°, but with the 
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same reservation as discussed in the previous section. Figure 17 shows 
the variation in dV/V (%) with change in V for two typical curves taken 
from Dodge (1934, plates IV and V). In addition, a portion of the frac- 
tional deviation curve for the “‘normal” case is included (upper right of 
Fig. 17). Where the uncertainties in ¢ are the same it is seen that greater 
accuracy can be obtained in the computation of V with small ¢ than is 
otherwise the case.!3 The curves D; and Dy» described in Fig. 16 are re- 
plotted in Fig. 17 for comparison of direct and computed values of V. 
As they cover about the same field in the diagram as the two curves 
@=10° and $=19°, it is seen that there is little to choose between the 
two methods under these conditions. If, however, the orientation condi- 
tions indicate a lower value of d¢ than that on which Fig. 17 is based, the 
computed value may be more accurate than that obtained by direct 
measurement. In all cases, direct measurement is preferable to the com- 
puted value based on = 45°. 

INTERFERENCE FIGURE METHODS. A number of methods based on in- 
terference figures are known (Johannsen 1918), but few are in wide- 
spread use at present. Where 2V<60° and a centered acute bisectrix 
figure is obtainable, Mallard’s method is convenient to use. The basic 
relation is sin V=D/K®§, where D =one-half the distance between the 
points of emergence of the optic axes of the crystal plate, K=a constant 
of the optical system (previously determined using minerals of known 
optic angle), 8=intermediate refractive index of the crystal under ex- 
amination. 

Assuming any random error in K to be negligible relative to the un- 
certainty in measuring 2D, differentiation of the above equation gives 

dV (radians) = tan V (> - =), 
Dewi 

If K=30 for the given optical system the uncertainty in 2D need not 
exceed 0.5 scale divisions, so that dD =.25. If d@ is known to +.001, it 
is easily shown that d@/@ is negligible relative to dD/D, so that, finally, 
dV (degrees) =57.3dD/D tan V. Curve M in Fig. 18 is based on this 
equation. It is seen that the accuracy of the method is low for small optic 
angles. 

Curve C shows the variation of dV/V with V where the optic angle is 
estimated from curvature of a single isogyre. With practice the uncer- 
tainty in estimating 2V by this method" need not exceed 6°, so that 
dV = 3°. Curve C is based on this value. 


18 The same comment made in the previous section regarding uncertainty in x for low } 
values applies here. 


M4 FE. S. Larsen, personal communication. 
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For comparison, curve D, (showing the fractional deviation in the di- 
rect method for dV= 1°) is included in Fig. 18. 

COMPARATIVE ACCURACY. Inspection of the diagrams (Figs. 14-18) 
makes apparent the difficulty of general comparison of optic angle meth- 
ods. Except where V and B are small, the refractive index method is 


Stine qjaresaG0e 45° 


Fic. 18. Interference figure methods. Variation in accuracy of V for Mallard’s method 
(curve M), where K=30 and dD=.25. Curve C represents the estimation method where 
dV =3°. Curve D, is for the direct method where dV=1°. 


undoubtedly the most accurate if the indices are determined by the 
minimum deviation prism method. However, where the indices are de- 
termined by matching with an oil to +.001, the method yields results 
which are not much better on the average than results obtained from 
estimating the curvature of an isogyre. 

Where favorable conditions for rotation exist, and with the uncer- 
tainty in ¢ not in excess of 1°, the extinction angle method will give 
values of V almost as free from error as the refractive index method (for 
dn=.0001 and B>.02). Where the uncertainty is +2°, however, the 
accuracy is approximately the same as the refractive index method with 
dn=.001. 

Comparison of the compensator method with other methods requires 
the assumption of a definite thickness for the section. Thus for t=.02 
mm. the dV/V compensator values are about intermediate between those 
obtained from the extinction angle and refractive index methods. 
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Although of restricted application, Mallard’s interference figure 
method is more accurate than direct measurement of optic angle with 
the universal stage, if the maximum values for dD and dV used in the 
diagrams are accepted as probable. 

Comparison of the three indirect universal stage methods with the 
direct method shows in most cases that the latter is at least equal, or 
definitely superior, in accuracy to the former. Exceptions occur with the 
refractive index method where dn=.0001 and B and V are moderate to 
large, and with the extinction angle method where ¢ is small and x large. 

For all methods less accuracy is possible for low values of optic angle 
than for high values. 


SUMMARY 


The usefulness of the diagrams presented in this paper depends on 
good judgment in selecting a reasonable maximum value of the un- 
certainty in any given measurement. In choosing the values on which the 
preceding figures are based the authors have assumed a skilled operator 
who has at his disposal good material and equipment. Where a series of 
readings of the same measurement is taken and a “‘probable error” com- 
puted, greater accuracy should be attained than shown here. Where the 
working material is poor less accuracy is to be expected. In general, 
where the uncertainties in measurement are believed to differ from those 
used in the diagrams, it is necessary only to substitute in the differential 
equation the revised value for the original and to construct a new graph. 

A set of these graphs readily available in the laboratory should prove 
highly useful by providing instant information regarding the order of 
magnitude of the maximum error to be expected in a given procedure. It 
is the authors’ hope that in time a more realistic attitude will be devel- 
oped toward mineralogic measurements than is now apparent from much 
of the published quantitative data. 
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DIAGNOSTIC CRITERIA FOR CLAY MINERALS* 
W. F. Brapiey, State Geological Survey, Urbana, Illinois. 


ABSTRACT 


The value of the utilization of complex formation between clay minerals and suitable 
types of organic liquids in the identification and characterization of clay minerals is illus- 
trated by means of a series of such complexes between ethylene glycol and various clay 
minerals. Particular attention is devoted to mixed layer minerals involving intergrowth of 
montmorillonite with illite. 


Recent observations on molecular associations between organic liquids 
and the montmorillonite type of clay minerals (1) have suggested that 
similar associations may be used to classify the other clay minerals. 
Not only is it true that many organic materials are available which as- 
sociate with montmorillonites in a clearly characteristic manner, but 
there is an extensive list of incompletely characterized clays in which low- 
temperature water seems sufficiently important that one might anticipate 
similar associations. 

For the limited number of observations reported here, the following 
thoughts have guided the choice of illustrative examples. The associa- 
tions of poly oxygen compounds with montmorillonite are assumed with- 
out the complicating factor of base exchange, and of these, the complex 
that contains ethylene glycol, the simplest molecule, invariably gave an 
excellent clearly characteristic x-ray diffraction diagram. Ethylene gly- 
col is both small enough to minimize possibilities of steric hindrance 
and sufficiently non-volatile to keep clay specimens damp during the 
registration of diffraction patterns, even though no complex may have 
been formed. 

THE MONTMORILLONITE GROUP 


Associations of montmorillonites with organic materials have been 
studied from time to time (2, 3, 4, 5) and have been the basis of Hen- 
dricks and Alexander’s color reaction for the group (6). X-ray diffraction 
diagrams represent only one of the methods by which montmorillonites 
are habitually recognized. Justification for the further complication of 
utilizing x-ray diffraction from chemically treated preparations must 
come from the characteristically more clean-cut appearance of patterns. 
Natural minerals of the montmorillonite group afford diffraction patterns 
which vary enormously in the relative prominence of the basal reflections 
(commonly used as the diagnostic criterion) and of the prismatic reflec- 
tions, to say nothing of the variations in diffraction angle for the first 
order of the base, and the vast differences in the perfection of develop- 
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ment of succeeding orders. Montmorillonites soaked in glycol and drained 
dry before registration of the diffraction diagram have invariably af- 
forded a clear sequence of higher orders of diffraction from the base, of 
greater intensity relative to prismatic reflections than are observed for 
natural materials, and representing a c-axis length not appreciably af- 
fected by previous history or humidity. This configuration has been in- 
terpreted to represent individual montmorillonite layers interleaved with 
double layers of glycol molecules, comprising a periodicity along cof 


nites compared with 
bentonite; (b) pink montmorillonite 
Mexico, U.S. Nat. Mus. no. 101836; (c) montmorillonite from 
and treated material below for each. 


Fic. 1. X-ray powder diffraction diagrams of natural montmorillo 
the diagrams after glycol treatment. (a) Wyoming 
from Tatatila, Vera Cruz, 
Rodeout, Utah. Natural clay above, 
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17 A. Diagrams have looked about equally clear for glycol complexes 
prepared from natural materials for which, when untreated, the higher 
orders of the base were: (a) clear and related to a 12.5 A periodicity; (b) 
clear and related to a 15 A periodicity; or (c) poorly defined or missing. 
Illustrative examples are shown in Fig. 1. In cases where montmoril- 
lonite is only a minor component in the presence of other material, dif- 
fraction effects from the complex are also more easily discerned than are 
those of the untreated material. 


THE KAOLINITE GROUP 


The kaolinite group of minerals, kaolinite, nacrite, and dickite, in- 
clude no molecular water in their crystallizations, and exhibit no re- 
actions with liquids which can be followed by diffraction methods. It 
is undoubtedly true that there is strong adsorption of organic liquids on 
the crystallite surfaces, as shown in the properties of materials like 
plasticene (or modelling clay), but the actual crystallization (and dif- 
fracticn patterns) are unaltered. 


ENDELLITE AND HALLOYSITE 


An unfortunately confused state of nomenclature with regard to hal- 
loysite has recently been resolved by the proposal of ‘‘endellite” as a 
name for the hydrated member (7). 

Both minerals are recognized readily enough on the basis of their nor- 
mal diffraction diagrams, but their behavior toward the organic liquids 
is interesting. The displacement of molecular water from endellite by 
ethylene glycol (and by other solvents) both serves to confirm its struc- 
ture and is an unmistakable criterion of identity. 

The characteristic strong diffraction line of endellite at 10.1 A is dis- 
placed to 10.8 A in the ethylene glycol complex. First, third, and fifth 
orders from this spacing, which is the complex thickness, are apparent 
in each case. These thicknesses represent the thickness of a kaolinite 
layer plus, respectively, one water layer of about 3.0 A and one glycol 
layer of about 3.6 A. In contrast with the montmorillonite complexes, 
which have two layers of glycol, it is pointed out that endellite has only 
one oxygen surface and montmorillonite has two oxygen surfaces availa- 
ble. It is assumed that the stability of each configuration is governed by 
the attraction of the oxygen surfaces for the slightly active methylene 
groups. 

The dependence of the mean refractive index of endellite on the con- 
stitution of the index oils was noted by Correns and Mehmel (8) but was 
questioned by Faust et al. (7). It is clear that the indices of complexes of 
the sort here described must differ appreciably from those of the natural 
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mineral. The calculated mean index for endellite is 1.49, and an esti- 
mated index for the glycol complex, based on the molecular volume of 
liquid glycol is 1.52, probably + 0.01. The observed index for the complex 
is 1.525, determined by rapid comparison with liquids made up from 
mineral oil and monochloronaphthalene. The endellite for this prepara- 
tion was selected from material from Eureka, Utah, and did not contain 
more than 5 or 10 per cent of halloysite. Complexes of endellite with 
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Frc. 2. The association of ethylene glycol] with halloysite. (a) natural endellite; (b) 
endellite treated with ethylene glycol; (c) halloysite resulting from aging of (a); (d) halloy- 
site treated with ethylene glycol; and (e) a plot of Hendricks and Teller’s function for 
units of 10.8 and 7.2 A, plotted to the approximate scale of the diffraction diagrams. The 
arrow indicates the position of normal reflection. 
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other active liquids (1) would also have indices ranging from 1.52 to 1.55. 

Halloysite derives from endellite upon loss of the molecular water. The 
resultant material is very porous, retains its original bulk, and readily 
absorbs water, but the structural change from the endellite crystalliza- 
tion to a kaolinitic crystallization is irreversible, or at best incompletely 
reversible. Wet halloysite specimens give essentially the halloysite dif- 
fraction diagram. 

The diffraction diagram of halloysite wet with ethylene glycol differs 
in a distinctive way from those of endellite, halloysite, or kaolinite. The 
prismatic interferences of endellite and halloysite are unchanged, but of 
the possible interferences related to the base, only one line, at 3.6 A, ap- 
pears sharp and well-defined. In the region corresponding to periodicities 
from about 11 to about 7 A. there appears only a broad weak band, in 
which the intensity seems to be somewhat higher near its edges. 

The function derived by Hendricks and Teller (9) for mixed layer 
minerals applies only to equally probable units of the same form factor. 
In the absence of suitable information about the form factors it is im- 
practical to attempt to compute a theoretical intensity distribution for 
this case (as will later be done for a more apt example), but a comparison 
of the plot of Hendricks and Teller’s function with the diffraction dia- 
gram in Fig. 2 indicates clearly that a random or near random interleav- 
ing of 10.8 and 7.2 A units would produce the same sort of diagram as is 
observed. Actually, it seems indicated that simple 7.2 A units outnumber 
the complex 10.8 A units, but the characterization of halloysite as the 
collapsed structure subject to random penetration of active liquids, seems 
justified. 


ILLITE 


Although illites invariably contain more or less molecular water, as 
is apparent from thermal analysis curves (10), the crystallization is de- 
termined by mica principles, and the position in the structure of any 
water not merely absorbed on exterior surfaces is not known. Treatments 
of illite with organic liquids do not therefore result in any structural 
changes, and diffraction diagrams do not differ from those of the un- 
treated materials. The frequently encountered instances in which illites 
are observed to be mixed with extraneous layers are discussed in connec- 
tion with the mixed layer minerals. 


Mixep LAYER MINERALS 


Several instances have been reported (11, 12, 13, 14, 15) in which it 
seemed that a specimen under consideration was actually composed of 
irregularly alternating layers of two distinct minerals. The treatment of 


DIAGNOSTIC CRITERIA FOR CLAY MINERALS 709 


diffraction effects from such mixed layers by Hendricks and Teller (9) 
has afforded a basis on which one can hope to examine in detail some of 
these suspected instances. 

A particularly fortunate case has been that of some bravaisite from the 
type locality, U. S. Nat. Museum Specimen No. 4918.* A study has 
shown not only that this material is about equally divided between two 
species, but also that the two form factors are similar enough to permit 
a reasonably accurate construction of the calculated diffraction effects. 

Grim and Rowland (10) concluded from thermal analysis curves that 
this material included both illite and montmorillonite. If it be assumed, 
then, that this material be made up of a random stacking of equaliy prob- 
able illite and montmorillonite layers, one case can be developed for the 
natural material itself, and a second case can be developed for the mix- 
ture of illite with glycol-complexed montmorillonite which would result 
from glycol treatment. 

In connection with the previous study of complexes of montmoril- 
lonite with various types of organic liquids (1), observed 00/ diffraction 
sequences from a series of complexes of varying cell heights have been 
available. Estimated intensities for the powder diffraction rings for this 
entire group have been plotted together against sin 6/A, and Fig. 3 is a 
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Fic. 3. Compositite curve of estimated relative intensities of 00/ powder diffraction 
rings versus sin #/A taken from a group of organic molecule-montmorillonite complexes of 
various c-axis periodicities. 


smooth average curve summarizing the composite data. Although each 
individual diagram has been rather strongly influenced by the organic 
matter, the allowable degree of error for the present purposes is so great 
that this curve adequately represents not only the contribution of the 
montmorillonite layers, but also the contribution of the similarly con- 
stituted, but rigid, illite layers. 

To arrive at the angular distribution of scattered radiation from the 
random mixed layer sequences assumed, then, it is only necessary to 


* Furnished through the courtesy of Dr. W. F. Foshag. 
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multiply the ordinates of Hendricks and Teller’s function for each case 
by the ordinates of Fig. 3. In Fig. 4 curves derived in this way for a mix- 
ture of 10 A and 15 A units, and of 10 A and 17 A units are compared 
with the diffraction diagrams of bravaisite and of glycol treated bravais- 
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Fic. 4. Bravaisite, Noyant Allier, France, U. S. Nat. Mus. no. 4918 (above) and the 
same treated with ethylene glycol (below), each compared with its synthesized curve for 
random mixed layers to scale approximately equivalent to the diffraction diagram. Erect 
arrows indicate positions of normal or near normal reflection. Inverted arrows indicate nor- 
mal reinforcements in the mixing function which are extinguished by the form factor. 
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ite, respectively. Reasonable agreement in the position, relative inten- 
sity, and apparent breadth of all the important maxima is evident. 

No other instance has been encountered in which the relative frequen- 
cies of illite and montmorillonite layers are nearly enough equal to con- 
form to Hendricks and Teller’s equation. However, the principle can 
probably be extended in a qualitative way to cases in which one species 
predominates over the other to a moderate extent. If for example, illite 
occurred more frequently than montmorillonite, we would expect from 


Fic. 5. Powder diffraction diagram of ‘‘glimmerton,” Sarospatak, Hungary, from U. 
Hofmann, treated with ethylene glycol (above), compared with illite purified from shale 
from Alexander County, Illinois, after the same treatment (below). 


inspection of the function that in general maxima would approach the 
position of the sharp interference maxima of illite. For the natural mate- 
rial, where the form factor extinguishes one of the small-angle maxima of 
the complex function, the diagram would only more nearly resemble illite. 
For glycol-treated material, however, both of the small-angle maxima 
should persist and should merely approach each other. 

Among the natural materials in which two small-angle maxima are 
developed by glycol treatment are several which have heretofore been 
described as illites or at least as mica-like clays. These include the “glim- 
merton” from Sarospatak, Hungary (16), the illite purified from an un- 
derclay in Grundy County, Illinois (17), a “metabentonite” from High 
Bridge, Kentucky (10), and possibly some glauconites. In Fig. 5 the dia- 
gram for glycol-treated “‘glimmerton” is reproduced. For specimens with 
lesser proportion of montmorillonite layers the diagrams are less clear. 
The ratio of mica to montmorillonite layers in the “glimmerton” is prob- 
ably near 2; in the other cases, somewhat higher. A diagram of a glycol- 
treated illite, unmixed with extraneous layers, is included for compari- 


son. 
Mixed layer minerals of other types also no doubt exist, of which par- 
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ticularly a mixture of montmorillonite and illite with montmorillonite 
predominating should be amenable to interpretation by this method, 
but no clear cut occurrence has yet been recognized. 


MISCELLANEOUS MINERALS 


The number of hydrous clayey or clay-like minerals which have been 
described and named is too great to attempt to catalogue them here. 
Probably the greatest number of these are the varieties which are actu- 
ally members of the montmorillonite group. Most of the specimens of 
this sort reported on by Grim and Rowland (10) have been made availa- 
ble for examination in this study. The characteristic 17 A glycol complex 
is readily prepared from all of the so-called beidellites discussed by them, 
except their number 8 A (U. S. Nat. Museum no. R7595), from the non- 
tronites, the white Hector magnesium clay, smectite, chloropal, and 
volkonskoite. 

Loosely related to these materials are two other mineral types in 
which comparable structural units appear, the vermiculites and atta- 
pulgite. 

Vermiculites—Although the ideal structure of vermiculite (12) im- 
mediately suggests that the interlayer water could he displaced by or- 
ganic liquids, no attempts to make such complexes from the familiar 
crystalline materials have yet been successful. However, a series of 
cryptocrystalline weathered amygdules collected by Dr. R. M. Grogan 
from basic Keweenawan flows along the north shore of Lake Superior, 
which have been identified as vermiculites by «-ray diffraction methods, 
do form the same sort of 17 A glycol complex as is observed for montmoril- 
lonites. In view of the known tendency for vermiculites to occur mixed 
with chlorites, it seems quite possible that such mixing, particularly if it 
were lateral within layers, might inhibit the introduction of glycol into 
the natural large crystals even though it would be readily introduced into 
a fine-grained, or into an ideal specimen. 

Attapulgite.—The association ofany individual water molecule with the 
silicate portion of attapulgite is probably due to hydrogen bonding, 
analogous with that exhibited by the water associated with montmoril- 
lonite, but the disposition of the silicate chains (18) confines the water 
to narrow channels and precludes the changes in dimensions that are 
observed for swelling clays. The water channels have cross sections of 
about 3.76.0 A, which might accommodate single strings of ethylene 
gycol, but in the absence of dimensional changes, only minor intensity 
variations are to be expected upon glycol treatment, and no obvious 
variations were noted in the diffraction diagrams. Refractive indices, 
however, are sharply reduced. Well-oriented flakes with y=1.540, 
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y—a=.032, showed after treatment y=1.50 with comparable bire- 
fringence. The calculated mean index for the silicate skeleton above is 
about 1.44. The observed indices can therefore apparently be considered 
to arise from extensive extraction of the hydrogen-bonded water accom- 
panied by only very limited penetration of glycol into the channels. 


SUMMARY 


The association of certain organic liquids with clay minerals affords 
a basis for the classification of such minerals into the principal groups, 
and at the same time serves to confirm the general features of the respec- 
tive structures. 
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ADDITIONAL OPTICAL AND CHEMICAL DATA ON THE 
STILPNOMELANE GROUP OF MINERALS 


C. OsBoRNE HuTTON, Geological Survey, Wellington, New Zealand. 


SUMMARY 


Since publication of a paper on the stilpnomelane group of minerals (Hutton, 1938), 
the author has succeeded in obtaining a member of this series with a composition consider- 
ably different from those previously analysed; in addition two analyses of stilpnomelanes 
have been published by other authors (Wager and Deer, 1939, p. 188; Ayres, 1940, p. 432), 
together with their optical constants. These analyses are discussed and the extent to which 
the provisional curves published earlier (Hutton, 1938, p. 187) require modification, is in- 
dicated. 


OpTICAL DATA 


The stilpnomelane newly analysed by Mr. F. T. Seelye for the writer 
was separated from an albite-epidote-stilpnomelane-chlorite schist con- 
taining minor amounts of quartz, titanite, muscovite, apatite, and cal- 
cite. The rock is coarsely crystalloblastic, strongly laminated, and 
schistose; it is a typical member of the green-schist facies as developed in 
the Chl. 4 subzone of the chlorite zone of western Otago (Hutton and 
Turner, 1936; Hutton, 1940). The stilpnomelane itself, in very coarsely 
crystalloblastic plates, and radiating or sheaf-like aggregates up to 5 mm. 
in length, is intimately interlaminated with plentiful chlorite from which 
it seems to be in the process of developing; it appears to lack any marked 
preferred orientation, and the sheaf-like aggregates in some cases lie 
transverse, in others parallel, to foliation and schistosity planes. 

Optical properties of the associated chlorite are as follows: 

a=1.635+0.002 
B=7=1.640 
y—a=0.005 
2V=0° 
Sign: negative 
Elongation: positive 

a=very pale greenish yellow to colourless. 

B=~y=deep grass green. 
Absorption y=6>a 
Dispersion r>v (not strong). 


Pleochroism 


These data indicate a composition just inside the field of brunsvigite 
in the sense used by Winchell (1936, p. 649, Fig. 4) in his most recent 
study of the chlorite group. In an earlier investigation of the chlorites 
occurring in schists of the chlorite zone of western Otago, the @ values of 
over thirty chlorites were found to be within the limits 1.607-1.639 (Hut- 
ton, 1940, p. 19). The chlorite just described is thus a relatively iron-rich 
type. 
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The stilpnomelane mineral was separated from the schist by careful 
crushing, screening, and centrifuging methods following the procedure 
recommended previously (Hutton, 1938, p. 183). The analysis under- 
taken by Mr. F. T. Seelye is set out in Table 1A. 

When a and y for the three newly analysed stilpnomelanes (Table 1, 
A, RB, C) are plotted against the ratio (Fe, Mg, Mn) O/ (Fe, Al)2O3 on 
the original diagram (Hutton, 1938, p. 187), the values given for the 
Otago mineral and that described by Ayres and Park (analyses A and B) 
agree reasonably well with the existing curves. In the case of the Green- 
land stilpnomelane (analysis C), however, the values for a and y¥ cer- 
tainly fall close to the extension of these curves, as stated by Wager and 
Deer (1939, p. 190), but only if the composition is not taken into account. 
However, in the case of the Greenland stilpnomelane, the values of a and 
y quoted indicate a (Fe, Mg, Mn) O/ (Fe, Al).O3 ratio (molecular) of 
about 30/70, whereas the ratio calculated from the corresponding analy- 
sis is approximately 52/48. This discrepancy is possibly due to difference 
in composition between the unanalysed material used for refractive index 
determinations, and that selected by Wager and Deer for analysis. Fur- 
thermore, the relatively high CaO content of the analysis (4.28 per cent— 
a figure much higher than that in any other analysis available to the 
writer) points to possible impurity of the analysed material. 

The writer now suggests a slight modification of the previously pub- 
lished curves to take into account the new data given in Table 1, A and 
B. This has been done in Fig. 1. Points for y of the refractive index el- 
lipsoid tend to lie along a straight line with no one point far removed 
from it, whereas with the exception of Nos. 2 and 11, none of the a values 
are far removed from the line representing a of the refractive index el- 
lipsoid. For stilpnomelane minerals either 6 or y can be determined with 
reasonable accuracy so long as thin, fully transparent platelets are used. 
Greater difficulty attends the determination of a owing to the micaceous 
habit of the mineral, and considerable care is required in order to obtain 
the minimum value. For the stilpnomelane mineral represented by analy- 
sis No. 4, the value for a is given by Gruner (1937, p. 914) as 1.58+0.01; 
this limit of accuracy is unusually wide, so that the possible range has 
be indicated in Fig. 1. 


CHEMICAL DATA 


No further information of the chemical properties of members of the 
stilpnomelane group is available. However, Ayres’ (1940, p. 433) state- 
ment that stilpnomelane from Crystal Falls was “not at all attacked by 
sulphuric” acid is surprising in view of the results of earlier experiments 
by the present writer (Hutton, 1938, p. 190). The newly analysed Otago 
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Numbers in parentheses refer to the analyses plotted in Fig. 5 of Hutton (1938, p. 187). 
1 (15). Western Otago. C. O. Hutton (1938, p. 184, anal. C, Table 3). 
2 (1). Mesabi Range, Minnesota. F. F. Grout and G. A. Thiel (1924, p. 230, anal. No. 

1). 

3 (8). Theodor mine, Aumenau, Lahn. J. Holzner (1933, p. 216, anal. No. 4). 
4 (12). Baern, Moravia. J. W. Gruner (1937, p. 913, anal. No. 4). 
ys Western Otago. C. O. Hutton (1938, p. 184, anal. D, Table 3). 
6 (4). Pen-y-rallt, Merionethshire. Hallimond (1924, p. 194, anal. No. 1). 
7 (10). Zuckmantel, Austrian Silesia. C. O. Hutton (1938, p. 184, anal. E, Table 3). 
8 (17). Anna Mine, Baern, Moravia. C. O. Hutton (1938, p. 184, anal. F, Table 3). 
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stilpnomelane (A in Table 1), when treated for 6 hours with 5N H.SO, 
on the water-bath was completely decolourized, leaving colourless, iso- 
tropic plates of hydrated silica, with a refractive index of 1.438 after 
drying at 100°C. for 15 minutes. 


TABLE 1. ANALYSES OF STILPNOMELANE 


A. B. Cc 

SiO: 44.52 42.42 45.29 
Al,O3 7.19 6.71 Soy! 
Fe.0; 27.32 33.24 23.95 
FeO SJoSiil 0.85 8.99 
TiO: 0.10 — 0.02 
MnO 0.42 2.27 1.14 
MgO 5.63 5.20 3.30 
CaO 0.23 _— 4.28 
Na2O 0.32 — nt.dt. 
K20 ilar) = nt.dt 
V203 nt.fd. —_— — 
Cr2O3 nt.fd. —< _ 
H20+ 6.70 8.33 6.12 
H.0— 2.70 1.45 1.79 

100.21 100.47 100.45 
a= 1.612+0.002 1.634 1.626 
B=y= 1.700 1.730 1.745 
y-a= 0.083 0.096 0.119 
2V 0° 0° 0° 
Pleochroism: 

a= Bright golden yellow As for 8 but weaker Golden yellow 
B=y= Very deep reddish Golden yellow to Dark reddish brown 

brown chestnut brown 
Absorption y=B>a y=B>a y=B>a 
Sp. Gr. 2.82 — = 
(Fe, Mg, Mn)O 44.47 38.70 51.88 
(Fe, Al)203 5) 68 61.30 48.12 


A. Stilpnomelane from albite-epidote-stilpnomelane-chlorite schist, p. 2264. Cowcliff Hill, 
Kawarau S.D. Analyst: F. T. Seelye. 

B. Stilpnomelane from ferruginous Upper Huronian slates, Crystal Falls, Michigan. 
Analysts: V. L. Ayres and B. Park. The following elements were found spectrographi- 
cally: Na, Li, Co, Ni, V, Ba, Zn. (V. L. Ayres, 1940). 

C. Stilpnomelane from spherical vesicles in granophyre inclusion. Mellemé, Kangerdlugs- 
suaq, E. Greenland. (L. R. Wager and W. A. Deer, 1939, pp. 188-189). 


9 (9). Western Otago. C. O. Hutton (1938, p. 184, anal. B, Table 3). 
10. Western Otago. Anal. A, Table 1, this paper. 
ie Crystal Falls, Michigan. V. L. Ayres (1940, p. 432); or anal. B, Table 1, this paper. 
12 (16). Western Otago. C. O. Hutton (1938, p. 184, anal. G, Table 3). 
The solid triangles represent the plot of typical parsettensite (J. Jakob, 1923). 
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NOTES AND NEWS 
PROBERTITE FROM LOS ANGELES COUNTY, CALIFORNIA 


JosepH Murpocn, University of California at Los Angeles. 


A small specimen sent the writer recently, marked “Lang Borax 
Mine,” proved on examination to be probertite. Since this mineral has 
been reported only from the Ryan and Kramer districts, in California, 
advantage was taken of a class field trip to the Lang area, where a good 
amount of the material was collected, and the occurrence confirmed. 
Identification was made by optical properties, the indices of refraction 
corresponding with those of probertite, and chemical tests confirmed 
these. 

The appearance of the mineral is like that described by Foshag! for 
ulexite from this locality. The mineral occurs in lenticular aggregates of 
groups composed of slender radiating prisms up to 15.5 mm. in dimen- 
sions, these groups themselves being more or less flattened. None of the 
specimens were found in place, as the material was all collected on the 
dump, but some with attached shale indicated their occurrence as lenticu- 
lar nodules in a shale matrix. 

This identity of appearance with Foshag’s material, and with a similar 
occurrence of ulexite described from the Suckow mine in the Kramer dis- 
trict,” led to the thought that the determination might have been in er- 
ror, and the mineral really ulexite. However, a careful re-check of the 
indices, and a determination of the water of crystallization confirmed its 
identity as probertite. Duplicate water determinations were run, giving 
values of 25.7% and 25.8% respectively, as compared with the theoreti- 
cal value of 25.63%. All doubt was definitely cleared up by the discovery, 
among the specimens, of a few which were clearly different from the oth- 
ers, and were actually ulexite. This ulexite is essentially similar in ap- 
pearance to the probertite, but is finer textured, much silkier in appear- 
ance, and lacking the definitely radial structure of the latter. Under the 
microscope it shows the normal ulexite indices, and the typical fine tex- 
ture, and twinning. In further confirmation, two specimens were found 
showing the two minerals together, with ulexite covering the probertite 
rosettes, and enclosing residual patches of the latter mineral. The es- 
sential structure of the probertite groups was preserved, although not 
the detail. Asa matter of fact, most of the material was unaltered prober- 
tite, with a very small number of the specimens showing any ulexite. 


1 Foshag, W. F., The origin of the colemanite deposits of California: Econ. Geol., 16, 204 


(1921). 
2 Murdoch, J., Crystallography of ulexite: Am. Mineral., 25, 755 (1940). 
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Since the collection was made at random, this indicates the probable ac- 
tual proportions of the two minerals in the deposit. 

This relationship of ulexite and probertite induced the writer to look 
again at some specimens from the Kramer district, collected as ulexite. 
These were found to have been labelled “‘ulexite with residual kramer- 
ite,” affording another instance of the same sequence. Foshag,’ in writing 
of the Ryan District, suggested that ulexite had been altered to prober- 


Ng Wj borax 


Fic. 1 


tite although he did not find the two together anywhere. The basis for 
this idea was that the denser mineral (probertite) would be more likely 
to form from ulexite under the influence of pressure generated in the fold- 
ing of the layers in this area. This does not appear to have been the case 
in the two cases which the writer has observed, since both at Lang and at 
Kramer the opposite sequence is clear. 

At Kramer there are two examples of this relationship. First, the ra- 
diating groups or clusters of probertite described by Schaller under the 
name kramerite, were found on one of the dumps by the writer as ulexite 
pseudomorphs, with occasional probertite cores. Second, specimens from 
the Suckow mine, just west of the other occurrence, show the curious 
fungus-like growths, and satiny cross-fiber veins mentioned from the 
kernite deposit,® in layers of clay and borax. Most of these are milky 
white, uniformly fine fibered ulexite, but some show residual patches of 
probertite, as noted casually at the time of collection. This can be identi- 


* Foshag, W. F., Probertite from Ryan, Inyo County, California: Am. Mineral., 16, 
340 (1931). 

4 Schaller, W. T., Borate minerals from the Kramer District, Mohave Desert: U.S.G.S. 
Prof. Paper 158, 139 (1930). 

5 Schaller, W. T., op. cit., p. 139. 
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fied easily by its glassier luster and coarser appearance, and some of the 
patches are large enough to provide material for optical determination. It 
was clearly the original mineral of these occurrences, but has been almost 
completely replaced. The colemanite produced by the further alteration 
of the ulexite, as observed by Schaller,® was not seen in any of the writer’s 
material. 

The relationship of probertite to ulexite in these specimens is illus- 
trated by the accompanying drawings (Fig. 1), which represent the oc- 
currence of the two minerals. There are three types of these, all from 
Kramer: (a) divergent, cone-shaped masses, usually in solid borax; (b) 
the fungus like forms, in layers of clay, or of clay and borax crystals; and 
(c) more or less hemispherical aggregates, occurring in narrow veins of 
secondary borax which cut the older ore. 

The writer wishes to express his thanks to Mr. E. H. Peebles of San 
Gabriel, California, who sent in the specimen which led to the discovery 
of this new occurrence of probertite. Mr. M. Vonsen, of Petaluma, Cali- 
fornia, has informed the writer that he had recognized this material from 
Lang as probertite, some t me ago. 


6 Schaller, W. T., op. cit., p. 138. 


Dr. Ralph E. Grim, petrographer and head of the Division of Petrography of the IIli- 
nois Geological Survey, has been made petrographer and principal geologist in charge of the 
Geological Resources Section. This section comprises the following divisions: coal, oil and 
gas, industrial minerals, clay resources and clay mineral technology, groundwater geology 
and geophysical exploration, areal and engineering geology, stratigraphy and paleontology, 
and subsurface geology. In the Geochemical Section of the Survey Dr. William F. Bradley 
has been made chemist and head of the Division of X-ray and Spectrography. 


BOOK REVIEWS 


THEORY OF X-RAY DIFFRACTION IN CRYSTALS by Wi1t1am H. ZACHARIASEN, 
New York, John Wiley & Sons, (1945), V & 255 pages, 30 figures. Price $4.C0. 


This excellent book, written by one of the leading workers in the field, contains a com- 
plete mathematical treatment of the phenomena encountered in «x-ray diffraction in crys- 
tals. 

Chapter I starts with the macroscopic concept of a crystal and of crystal faces, intro- 
duces the lattice concept, the reciprocal lattice, and Fourier development of crystal prop- 
erties. 

Chapter II contains an excellent and concise discussion of the symmetry properties of 
crystals. Throughout the chapter dyadic algebra is used rather than geometric considera- 
tions. The properties of dyadics required, as well as the group theory needed, are derived 
in Appendix A and B. The possible symmetry operations in crystals are discussed. The 
properties of the 32 point groups are derived and discussed. All possible space groups are 
listed, and their complete derivation is treated. 

Chapter III deals with x-ray diffraction in ideal crystals. The Laue and Bragg equations 
are deduced, and the construction of the diffracted vector discussed. The main part of the 
chapter is devoted to the calculation, from the atomic scattering power, of the intensity dis- 
tribution in the interference maxima. Formulae are developed for the percent reflection, 
integrated scattering power, and half width in a number of cases. The cases treated are 
those of small crystals, thin, thick, and intermediate crystal plates, with and without ab- 
sorption. 

The agreement of these formulae with experiment is in general poor. The theory pre- 
dicts smaller integrated intensity and narrower intensity distribution than those found 
experimentally at the Laue-Bragg spots. The reason for the difference is to te found in the 
crystal imperfections. 

These imperfections are taken up in Chapter IV, where the discrepancy between ex- 
periment and theory is resolved. The mosaic structure of crystals is shown to be able to 
account for the broader and more intense interference maxima which are observed. Another 
important disorder is the heat motion. The normal coordinates of the lattice are introduced, 
and their thermal vibrations are shown to give rise to diffuse scattering and a decrease of 
intensity of the interference maxima. Temperature motion taken alone, without the im- 
perfections due to mosaic structure, would increase the disagreement ketween theory and 
experiment. The last section discusses some common disorders found in real crystals. 

The book is clearly and concisely written throughout. All the material in it is com- 
pletely and logically derived from the fundamental concepts. The numerous formulae for 
intensity and half width ought to be very useful for a worker in the field. A considerable 
amount of the material of the book consists in original work of the author and cannot be 
found in any publication. This is especially true of the treatment of mosaic structure in 
Chapter IV. 

The book is extensively mathematical. Dr. Zachariasen ends his introduction with the 
sentence: ‘“The modest size of this book is not a fair measure of the time I have spent on 
it.” This is indeed to be believed. The reader trying to master the material of this book may 
find the same statement to apply to himself. 

Maria GOPPERT MAYER 
Dept. of Chemisiry 
Columbia University 
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PYRITES DEPOSITS OF MISSOURI by Oxtver R. Grawe has recently been issued 
as volume 30 of the Missouri Geological Survey and Water Resources. This report of 482 
pages describes the location, history, geologic conditions and production of 91 mires and 
prospects. Spectrographic analyses of the pyrites and its oxidation products have shown 
the presence of 23 elements, most of which are present in amounts less than 0.05%. The 
report also includes x-ray analyses, studies of polished sections and a description of 36 min- 
erals that have been reported as occurring in the pyrites and iron ore bearing sink structures 
of the northern Ozark plateau. 


REPORT OF INVESTIGATIONS NO. 50 OF THE SOUTH DAKOTA GEOLOGI- 
CAL SURVEY embraces a “‘PRELIMINARY REPORT ON THE MINERALOGY OF SOME PEG- 
MATITES NEAR CusTER” by D. Jerome Fisher. The main object is to suggest an economical 
plan for the recovery of the valuable minerals in the pegmatites of the area. This report 
covers the history, geology, detailed mineralogy and paragenesis of five pegmatites desig- 
nated by the following Lodes: Buster, Earl, Custer Mountain, Old Mike and Ross. Descrip- 
tions of the geology of three of the pegmatites of the Custer District appear in Report of 
Investigations No. 44 of the South Dakota Geological Survey. 


PHILADELPHIA MINERALOGICAL SOCIETY 
Academy of Natural Sciences of Philadelphia, May 5, 1945 


Dr. W. Hersey Thomas presided, with 59 present. Dr. Hugh E. McKinstry of the For- 
eign Economic Administration addressed the society on ‘‘Collecting minerals for the war.” 
Sources of mica, beryl, iceland-spar, vanadium and tungsten were described. 


Meeting of June 7, 1945 


Dr. Thomas was in the chair, with 54 present. Dr. William J. Kirkpatrick addressed 
the members on ‘‘Growth of Crystals.’”’ The historical introduction to the growth of euhe- 
dral crystals included Robert Boyle’s (1669) observation that the rate of growth influences 
crystal habit, Romé Delisle’s (1783) observation that sodium chloride crystallizes in octa- 
hedra from urine, and Le Blanc’s (1808) observation that alum is octahedral when grown 
in weak alkali solutions. Beudant (1818) published a paper of 100 pages on the causes that 
determine the variations of crystalline forms. 

Some factors that influence crystal habit are pH of a solution, presence of cosolutes, 
colloids or dispersing agents, rate of growth, temperature, pressure, and degree of agita- 
tion. A law of crystal growth may be stated: “Tf all of the faces grow at the same rate, 
the faces will disappear in the order of their distance from the center, d=1/V//2+2?+2, 
where h, k, and / are the Miller indices. 

Sodium chloride usually crystallizes in cubes; the presence of urea and other alkalies in 
the solution cause it to crystallize in striated octahedra, although no record could be found 
of anyone having obtained good reflections on the goniometer from octahedral faces. 

Some practical applications of crystal habit control were given: the dodecahedral face 
is the most effective plane in nickel for catalyzing hydrogenation reactions. Because they 
dissolve at a predetermined rate, many commercial chemicals such as photographer’s hypo 
and sodium chromate are supplied in uniform equidimensional crystals. The habit of cupric 
chloride is influenced by the presence of certain tissue extracts: a fact made use of in pathol- 


ogy. 
J. S. FRANKENFIELD, Secretary. 
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NEW YORK MINERALOGICAL CLUB, INC. 
Abstract of Meeting of May 16, 1945 


The president and treasurer gave their annual reports and a trip to the William Boyce 
Thompson collection was announced. Mr. Newton Marshall, Vice-president of the South 
American Gold and Platinum Co. described dredging operations for precious metals in 
the Choco region of South America, illustrating his lecture with a motion picture. 

ELIzABETH ARMSTRONG, Secretary. 


Abstract of Meeting of October 17, 1945 

With an attendance of about eighty, the New York Mineralogical Club opened its fall 
season with the usual meeting on summer collecting by members. In spite of the difficulties 
of travel members managed to visit New Mexico, Colorado, Pennsylvania, New York, New 
Hampshire, Maine, Vermont, Connecticut, New Jersey, and Massachusetts localities with 
the usual varying degrees of success. The most notable finds were 18 faceted bery] crystals, 
good New Jersey carnelian, Connecticut green gem tourmalines, and minerals from St. 
Peter’s Dome, Colorado. 

Mrs. E. J. Marcin, Secretary. 


NEW MINERAL NAMES 


Alumino-chrysotile 


D. P. SERDyUCHENKO: Alumino-chrysotile, a member of the isomorphous series serpen- 
tine-parakaolinite. Compt. rend. acad. sci. U.R.S.S., 46, 117-118 (1945). 

Flattened-columnar green chrysotile, with fibers 15-20 cm. in length, was found in the 
basin of the Markopidj River, North Caucasus. Analysis gave: SiQ2 42.54, Al,O3 5.68, 
Fe203 1.06, FeO 0.74, MgO 35.57, CaO 0.13, H2O— 0.38, H20+ 13.26; sum 99.36%. This 
corresponds to 87% chrysotile, 13% kaolinite. The mineral occurs at the contact of granite 
with serpentinites. 

Discussion: An unnecessary name for aluminian chrysotile. 

MICHAEL FLEISCHER. 


Unnamed 


ANDERS BystROM: Monoclinic magnetic pyrites. Arkiv. for Kemi, Mineral. Geol., 19B, 
No. 8, 1-8 (1945). 

X-ray study was made of 19 pyrrhotite samples from Swedish localities. Five were found 
to be hexagonal, with a=3.434 to 3.437+.002 A, c=5.714 to 5.726+.003 A. Three samples 
were analyzed for S giving 51.93, 53.03 and 53.20.atomic % S. Ten samples were found to 
be monoclinic, with a= 5.933 to 5.941 A, b=3.425 to 3.430, c=5.677 to 5.689 A, B=89.56° 
to 89.64°. Nine samples were analyzed for S, giving 51.74 to 53.01 atomic % S. Four sam- 
ples were found to be mixtures of the hexagonal and monoclinic forms. A sample of the 
monoclinic mineral and one of the mixtures were transformed into hexagonal phases when 
heated for 15 hours at 600°C. in evacuated quartz tubes. 
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